AFIT/GE/ENG/97D-12 


EFFECTS  OF  NEAR-FIELD  SCATTERERS 

ON 

SPACE-TIME  ADAPTIVE  PROCESSING 
THESIS 

Jonathan  W.  Fitton 
Captain,  USAF 

AFIT/GE/ENG/97D-12 


Approved  for  public  release;  distribution  unlimited 

UHC  QUAUTY  mBPSUTSD  3 

101  ' 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  Department  of  Defense  or  the  United  States  Government. 


AFIT/GE/ENG/97D-12 


EFFECTS  OF  NEAR-FIELD  SCATTERERS 

ON 

SPACE-TIME  ADAPTIVE  PROCESSING 

THESIS 

Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 
In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science  in  Electrical  Engineering 

Jonathan  W.  Fitton,  B.S.E.E. 

Captain,  USAF 

December,  1997 


Approved  for  public  release;  distribution  unlimited 


AFIT/GE/ENG/97D-12 


EFFECTS  OF  NEAR-FIELD  SCATTERERS 

ON 

SPACE-TIME  ADAPTIVE  PROCESSING 


Jonathan  W.  Fitton,  B.S.E.E. 
Captain,  USAF 


Approved: 


^/C  Af!  OCT  9 

Capt  William  Melvin,  Ph.  D.  ^  Date 

Committee  Member 


Acknowledgements 

I  would  like  to  thank  my  advisor,  Dr.  Byron  Welsh,  for  introducing  me  to  space-time 
adaptive  processing  and  for  his  assistance  in  preparing  this  thesis.  I  would  also  like 
to  thank  my  sponsor,  Capt.  Bill  Melvin  of  Rome  Labs,  for  suggesting  this  topic. 
I  would  like  to  thank  my  committee  members,  Maj.  Pete  Collins  and  Maj.  Mike 
Temple,  for  their  inputs  in  preparing  this  thesis.  And  finally,  I  would  like  to  thank 
Maj.  Scott  Berger  for  helping  me  find  the  sign  error  in  my  derivation  that  had  me 
going  in  circles  for  weeks. 

Jonathan  W.  Fitton 


111 


Table  of  Contents 

Page 

Acknowledgements .  iii 

List  of  Figures  .  vii 

List  of  Tables .  x 

Abstract .  xi 

I.  Introduction  .  1 

II.  Multidimensional  Signal  Processing .  3 

2.1  Introduction .  3 

2.2  Multidimensional  Signals .  3 

2.3  Beamforming .  5 

2.3.1  Weighted  Delay-and-Sum  Beamformer  ....  6 

2.3.2  Filter-and-Sum  Beamforming .  11 

2.3.3  Frequency-Domain  Beamforming .  13 

2.4  Discrete  Time  Beamforming .  14 

2.5  Total  Antenna  Pattern .  17 

2.6  Summary .  19 

III.  Space-Time  Adaptive  Processing .  20 

3.1  Introduction .  20 

3.2  Radar  System .  20 

3.3  Target .  23 

3.4  Noise  .  28 

3.5  Jamming .  29 

iv 


Page 

3.6  Clutter  .  31 

3.6.1  Effects  of  Platform  Motion  on  Clutter .  34 

3.6.2  Rank  of  the  Clutter  Covariance  Matrix  ....  36 

3.6.3  Intrinsic  Clutter  Motion .  40 

3.7  Fully  Adaptive  Space-Time  Adaptive  Processing  (STAP)  43 

3.8  Signal-to-Interference-plus-Noise  Ratio .  46 

3.9  Summary .  46 

IV.  Near-Field  Scatterers .  49 

4.1  Introduction .  49 

4.2  Radar  System .  52 

4.3  Scattered  Target  Signal .  54 

4.4  Scattered  Jammer  Signal .  57 

4.5  Scattered  Clutter  Signal .  59 

4.6  Fully  Adaptive  Processing  in  the  Presence  of  Near-field 

Scatterers .  61 

4.7  Summary .  64 

V.  Results  .  69 

5.1  Test  Plan .  69 

5.2  Results  .  69 

5.2.1  No  External  Interference .  69 

5.2.2  External  Interference  Present  .  86 

VI.  Conclusions .  101 

6.1  Conclusions .  101 

6.2  Suggestions  for  Future  Research .  102 

Appendix  A.  Notational  Conventions  and  Matrix  Operations  ....  103 


V 


Page 

Appendix  B.  Supplemental  Figures .  106 

Bibliography  .  114 

Vita .  115 


vi 


List  of  Figures 

Figure  Page 

1.  Loci  of  points  in  (k,  /) .  5 

2.  Bandpass  filter  in  wavenumber-frequency  space .  6 

3.  Block  diagram  for  the  weighted  delay-and-sum  beamformer.  .  6 

4.  Lowpass  Filtering  of  a  misaligned  signal .  9 

5.  Array  pattern  for  a  10  element  linear  array  with  no  time  delays.  10 

6.  Effect  of  grating  lobes  in  (k,  /)-space .  11 

7.  Bandpass  filter  in  wavenumber-frequency  space .  15 

8.  Quantized  Steering  Delays .  16 

9.  Effect  of  steering  delay  quantization  on  array  pattern .  17 

10.  Radar  system  and  antenna  array  geometries .  21 

11.  Radar  CPI  datacube .  22 

12.  Relationship  between  wavenumber  vector,  spatial  frequency,  wave¬ 
length,  and  element  spacing .  24 

13.  Clutter  Geometry .  31 

14.  Clutter  covariance  matrix  =  1),  showing  the  Toeplitz-block- 

Toeplitz  structure .  35 

15.  Clutter  ridges  for  different  values  oi  (3 .  37 

16.  Clutter  ridges  for  different  misalignment  angles .  38 

17.  Clutter  ridges  for  different  combinations  of  (f>a  and  P .  38 

18.  Clutter  eigenspectra .  39 

19.  Effective  array  positions  for  successive  pulses .  40 

20.  Clutter  eigenspectra  different  misalignment  angles .  41 

21.  The  adapted  pattern  in  the  presence  of  interference .  45 

22.  Four  possible  scattering  paths  with  a  near-field  scatterer  ...  50 

23.  Geometry  of  the  radar  system  including  the  near-field  scatterer.  51 

vii 


Figure  Page 

24.  Location  of  the  near-field  scatterers  relative  to  the  array.  ...  61 

25.  Clutter  covariance  matrix  in  the  presence  of  seven  near-field 

scatterers .  62 

26.  Upper  left  corner  of  the  clutter  covariance  matrix .  63 

27.  The  adapted  pattern  in  the  presence  of  interference  and  near¬ 
field  scatterers . 65 

28.  Near-field  scatterer  orientations .  70 

29.  SNR  for  5  NFS,  orientation  A  .  71 

30.  SNR  for  7  NFS,  orientation  A  .  72 

31.  SNR  for  8  NFS,  orientation  A  .  73 

32.  SNR  for  15  NFS,  orientation  A .  74 

33.  Effect  of  eight  NFS  on  the  direct-path  wavefront  perceived  by 

a  fifteen  element  array.  .  75 

34.  Signal-to-noise  ratio  for  orientation  A,  no  interference  ....  77 

35.  Signal-to-noise  ratio  for  orientation  B,  no  interference  ....  78 

36.  Signal-to-noise  ratio  for  orientation  C,  no  interference  ....  79 

37.  Signal-to-noise  ratio  for  orientation  D,  no  interference  ....  80 

38.  Signal-to-noise  ratio  for  orientation  E,  no  interference  ....  81 

39.  Signal-to-noise  ratio  for  orientation  F,  no  interference  ....  82 

40.  Adapted  antenna  pattern  as  a  function  of  azimuth  angle,  ^^  =  0°.  83 

41.  Adapted  antenna  pattern  as  a  function  of  azimuth  angle,  (j)t  = 

20° .  84 

42.  Adapted  antenna  pattern  as  a  function  of  azimuth  angle,  (j)t  = 

40° .  85 

43.  Eigenvalues  of  the  clutter  covariance  matrix  for  different  orien¬ 
tations .  87 

44.  Adapted  antenna  pattern  for  orientation  A,  =  0° .  89 

45.  Adapted  antenna  pattern  for  orientation  B,  0*  =  0° .  90 

46.  Adapted  antenna  pattern  for  orientation  C,  (^t  =  0° .  91 


vm 


Figure  Page 

47.  Adapted  antenna  pattern  for  orientation  D,  =  0° .  92 

48.  Adapted  antenna  pattern  for  orientation  E,  (j)t  =  0° .  93 

49.  Adapted  antenna  pattern  for  orientation  F,  =  0° .  94 

50.  SINK  for  orientation  A,  <j)t  =  20° .  95 

51.  SINR  for  orientation  B,  <l>t  =  20° .  96 

52.  SINR  for  orientation  C,  <f)t  =  20° .  97 

53.  SINR  for  orientation  D,  =  20° .  98 

54.  SINR  for  orientation  E,  <f)t  =  20° .  99 

55.  SINR  for  orientation  F,  =  20° .  100 

56.  Angle-Doppler  pattern  with  no  scatterers .  107 

57.  Angle-Doppler  pattern  for  orientation  A .  108 

58.  Angle-Doppler  pattern  for  orientation  B .  109 

59.  Angle-Doppler  pattern  for  orientation  C .  110 

60.  Angle-Doppler  pattern  for  orientation  D .  Ill 

61.  Angle-Doppler  pattern  for  orientation  E .  112 

62.  Angle-Doppler  pattern  for  orientation  F .  113 


IX 


List  of  Tables 


Table  Page 

1.  Symbols  for  Radar  System  Parameters .  27 

2.  Parameters  used  to  generate  Fig.  21 .  44 

3.  Symbols  for  near-field  scatterer .  53 

4.  Definition  of  orientations .  69 

5.  Rank  of  the  clutter  covariance  matrix  for  each  orientation  and 

number  of  scatterers  S .  86 


X 


AFIT/GE/ENG/97D-12 


Abstract 

Near-field  scatterers,  such  as  the  wing  of  an  airplane,  can  affect  the  ability  of 
the  processor  to  null  out  clutter  and  jammer  signals.  Target,  clutter  and  jammer 
signals  will  refiect  off  the  near-field  scatterers  into  the  array  and  appear  to  be  coming 
from  a  direction  different  from  their  true  sources.  In  many  analyses,  it  is  assumed 
that  the  signals  entering  the  array  are  plane  waves.  This  assumption  will  not  hold 
for  the  scattered  signals,  which  will  need  to  be  treated  as  spherical  waves,  further 
complicating  the  computations. 

This  thesis  develops  a  theoretical  model  similar  to  that  in  [8],  but  including 
the  effects  of  near-field  scatterers.  The  optimum  weight  vector,  normally  computed 
using  the  steering  vector  to  the  target  and  the  covariance  matrix  of  the  undesired 
signals,  must  now  include  the  effects  of  the  scattered  signals  as  well.  This  thesis 
shows  that  the  space-time  steering  vector  for  the  scattered  signal  can  be  written  in 
a  form  similar  to  the  direct  path  signal.  The  total  space-time  snapshot  of  a  signal  is 
the  sum  of  the  direct  path  and  scattered  path  steering  vectors  associated  with  that 
signal. 

The  near-field  scatterers  cause  more  signal  energy  to  enter  the  array  than 
would  be  the  case  without  the  scatterers.  When  the  scatterers  are  more  than  a 
half- wavelength  apart,  they  create  grating  lobes  in  the  adapted  antenna  pattern  and 
a  loss  of  signal-to-interference-plus-noise  ratio  (SINK).  When  the  scatterers  are  less 
than  a  half  wavelength  apart,  they  tend  to  reduce  the  sidelobe  level  and  increase  the 
SINK.  The  rank  of  the  clutter  covariance  matrix,  which  is  an  indicator  of  the  number 
of  degrees  of  freedom  required  to  eliminate  the  clutter  return,  increases  dramatically 
when  near-field  scatterers  are  present. 
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EFFECTS  OF  NEAR-FIELD  SCATTERERS 

ON 

SPACE-TIME  ADAPTIVE  PROCESSING 


I.  Introduction 

Space-time  adaptive  processing  (STAP)  is  a  technique  to  enhance  the  target  detec¬ 
tion  capability  of  a  radar  array  by  adjusting  the  parameters  of  the  array  receivers 
to  produce  antenna  pattern  nulls  in  the  direction  of  undesired  signals,  such  as  clut¬ 
ter  and  jamming,  while  placing  the  main  lobe  of  the  pattern  in  the  direction  of 
the  target.  The  receiver  parameters  also  control  how  different  Doppler  frequencies 
are  processed.  Because  the  characteristics  of  the  target  and  undesired  signals  are 
unknown,  they  must  be  determined  from  the  received  signals.  Then  the  receiver 
parameters  are  adjusted  adaptively  to  form  the  desired  beam  pattern.  This  adaptive 
beamforming  in  the  spatial  and  temporal  domains  is  the  essence  of  STAP. 

The  fundamentals  of  adaptive  beamforming  were  first  published  by  Apple- 
baum  [1]  and  Widrow  [9].  In  1973,  Brennan  and  Reed  [5]  published  the  first  work 
describing  optimal  space-time  filtering  using  adaptive  arrays.  At  that  time,  the  com¬ 
putational  complexity  of  space-time  adaptive  processing  was  beyond  the  capability 
of  the  technology  that  existed  to  be  implemented  effectively.  Advances  in  digital  sig¬ 
nal  processing  since  that  time,  along  with  the  development  of  new  algorithms,  have 
made  implementation  of  STAP  in  real-time  applications  more  feasible.  In  1994, 
Ward  [8]  published  a  theoretical  development  of  space-time  adaptive  processing.  In 
his  report,  he  presents  models  for  each  of  the  signals  received  by  an  airborne  radar 
array  and  shows  how  they  are  processed  using  optimal,  fully  adaptive  processing. 
He  also  presents  several  suboptimal  processing  methods. 
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Most  works  published  on  STAP  assume  the  array  is  in  free  space.  They  do 
not  account  for  the  effect  the  airframe  carrying  the  array  has  on  the  received  signal. 
Specifically,  how  does  a  scatterer  in  the  near-field  of  the  array  affect  the  received  sig¬ 
nal  and  the  ability  of  the  processor  to  null  out  interference.  In  [4],  the  authors  show, 
through  simulation,  how  a  near-field  scatterer  can  decrease  target  detectibility  by 
causing  clutter  returns  at  the  same  Doppler  frequency  as  the  target  to  appear  to  be 
coming  from  the  same  direction  as  the  target.  Others,  such  as  Yu  et  al.  [10],  have  de¬ 
veloped  computer  simulations  combining  the  geometrical  theory  of  diffraction  (GTD) 
with  finite  element  modeling  (FEM)  to  compute  the  pattern  of  an  array  mounted  on 
a  specific  airframe.  However,  there  has  not  been  a  theoretical  development  similar  to 
Ward’s  which  accounts  for  a  near-field  scatterer,  and  this  development  is  the  purpose 
of  this  thesis. 

The  next  chapter  develops  the  fundamentals  of  beamforming,  following  the 
development  presented  in  [6].  Chapter  3  traces  through  Ward’s  development  of  fully 
adaptive  space-time  adaptive  processing.  Chapter  4  extends  Ward’s  derivation  to 
account  for  a  near-field  scatterer.  Chapter  5  presents  simulation  results  and  Chapter 
6  presents  the  conclusions  of  this  thesis.  Throughout  this  thesis,  matrix  notation  is 
used  to  simplify  expressions.  Appendix  A  describes  the  notation  and  conventions 
used  in  this  thesis. 
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II.  Multidimensional  Signal  Processing 

2.1  Introduction 

The  purpose  of  space-time  adaptive  processing  (STAP)  is  to  adjust  an  array  so 
that  an  incoming  desired  signal  can  be  detected  while  at  the  same  time  suppressing 
undesired  signals  from  jammers,  clutter  and  other  sources.  The  first  step  in  the  de¬ 
velopment  of  space-time  adaptive  processing  is  an  introduction  to  multidimensional 
signal  processing  and  array  beamforming.  This  chapter  introduces  the  concepts  of 
multidimensional  signal  processing  and  filtering  in  the  wavenumber-frequency  do¬ 
main.  It  then  develops  expressions  for  different  types  of  beamformers,  and  concludes 
with  a  discussion  of  the  effect  of  digital  implementation  of  these  beamformers. 

This  chapter  closely  follows  the  development  in  Chapter  6  of  [6]. 

2.2  Multidimensional  Signals 

A  space-time  signal  can  be  described  as  a  four-dimensional  (4-D)  function 
s(x,t),  where  x  is  a  three  dimensional  vector  indicating  position  and  t  is  time.  This 
function  can  be  analyzed  as  a  four-dimensional  Fourier  transform.  The  function 

5(k, /)=  /  /  s(x,t)exp[— j27r(/t  —  k^x)]dxdt  (1) 

J— 00  J  —  OO 

is  the  wavenumber-frequency  spectrum  of  the  signal.  The  variable  /  is  the  temporal 
frequency  of  the  signal,  and  k  is  the  wavenumber  vector.  The  wavenumber  vector 
represents  the  spatial  frequency  of  the  signal,  or  the  number  of  wavelengths  of  the 
signal  per  unit  distance  in  each  of  the  three  orthogonal  spatial  directions.  The  term 
k^x  is  the  inner  product  of  the  wavenumber  vector  k  =  {kx,  ky,  k^Y  the  position 
vector  X.  The  inverse  transformation  is  defined  as 

=  /  S{k,t)exp\j2n{ft-'klx.)]dkdf.  (2) 

)“  J—OO  J-oo 
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A  propagating  plane  wave  is  represented  by  an  elemental  signal  of  the  form 


e(x,  t)  =  exp{j2iT{fot  -  k^x)].  (3) 

An  alternate  way  of  writing  Eqn.  (3)  is 

e(x,  t)  =  exp[j27r/o(t  -  a^x)],  (4) 


where 


A  ko 


(5) 


From  this  definition,  it  can  be  seen  that  e(x,  t)  represents  a  plane  wave  propagating 
in  the  direction  ocq  with  a  speed  equal  to  1/1q!o|-  The  vector  olq  is  referred  to  as 
the  slowness  vector  because  its  magnitude  is  equal  to  the  inverse  of  the  propagation 
speed.  The  slowness  vector  is  a  convenient  way  to  describe  both  the  frequency  and 
direction  of  propagation  of  a  signal. 

Taking  the  Fourier  transform  of  (4)  produces 


B(k,/)  =  i(k-k„){(/-/„) 


(6) 


which  represents  a  four-dimensional  impulse  function  in  (k,  /)-space  at  the  point 
k=ko  and  /  =  /o-  Thus,  every  point  in  (k, /)-space  corresponds  to  a  plane  wave 
in  (x,t)-space  with  a  particular  orientation  and  frequency.  In  Fig.  1,  various  loci 
of  points  are  shown  for  three  dimensions  {kx,ky,f).  (The  k^  component  is  omitted 
to  simplify  the  figures.)  Figure  1(a)  shows  that  all  signals  with  the  same  frequency 
lie  in  a  plane  normal  to  the  /  axis.  All  signals  coming  from  the  same  direction  will 
lie  in  a  half-plane  normal  to  the  (kx,ky)  plane,  as  shown  in  Fig.  1(b).  Figure  1(c) 
represents  all  signals  with  the  same  propagation  speed.  These  signals  lie  on  a  cone 
given  by  |k|  =  f/c.  Generally,  it  will  be  assumed  the  speed  of  propagation  is  a 
constant  equal  to  c,  the  speed  of  light  in  free  space. 
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Figure  1  Loci  of  points  occupied  in  (k,  /)-space  by  plane  waves  with  the  same 
(a)  frequency  /,  (b)  direction  of  propagation  k  ,  and  (c)  velocity  of 
propagation  l/\otQ\. 

Given  a  signal  s(x,t),  the  response  f{x,t)  of  a  linear,  shift-invariant  filter 
h{x,  t)  is  given  by  the  4-D  convolution  integral 

/(x,0=/  /  h{x-^,t-T)s{tr),d^dT.  (7) 

00  J—oo 

Transforming  this  expression  to  the  wavenumber-frequency  domain  produces 


F(k,/)=S(k,/)S(k,/). 


(8) 


The  task  is  to  design  the  filter  If(k,  f)  such  that  it  passes  signals  from  the  desired 
regions  of  (k,  /)-space  while  eliminating  signals  from  other  regions. 


2.3  Beamforming 

Figure  2  shows  how  a  3-D  bandpass  filter  might  be  represented.  Signals  with 
frequencies  in  the  range  of  A/  and  directions  of  arrival  within  the  range  of  angles 
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Figure  2  Bandpass  filter  in  wavenumber-frequency  space. 


Figure  3  Block  diagram  for  the  weighted  delay-and-sum  beamformer. 

A7  would  be  allowed  to  pass  through  the  filter;  all  other  signals  would  be  filtered 
out.  The  process  by  which  this  is  accomplished  is  called  beamforming.  To  achieve 
this  effect,  an  array  of  N  elements  is  arranged  in  space  such  that  the  nth  element 
is  located  at  x„(n  =  0, 1, . . . ,  iV  —  1).  Since  the  relative  position  of  the  elements  is 
fixed,  they  sample  the  signal  s(x,t)  spatially.  The  signal  from  the  nth  receiver  will 
be  given  by  r„(t)  =  s(x„,t),  and  these  signals  can  be  processed  by  various  methods 
in  either  the  time  or  frequency  domain. 


2.3.1  Weighted  Delay-and-Sum  Beamformer.  The  weighted  delay-and-sum 
beamformer  of  Fig.  3  is  one  of  the  simplest  beamforming  systems.  The  output  bf(t) 
is  the  average  of  the  weighted  and  delayed  receiver  signals 


(9) 


bfit)  =  4  E  Wnrnit  -  Tn), 

n=0 

where  r„(i)  is  the  signal  received  by  the  nth  receiver  Wn  and  are  respectively  the 
weight  and  relative  delay  of  the  nth  receiver.  Adjusting  the  delays  r„  has  the  effect 
of  “steering”  the  center  of  the  beamformer’s  passband  to  the  desired  orientation  in 
(k,  /)-space.  The  receiver  outputs  due  to  a  signal  from  the  desired  direction  will 
added  together  in  phase,  while  receiver  outputs  due  to  signals  from  other  directions 
will  be  out  of  phase  and  will  tend  to  cancel  each  other  out.  To  pass  plane  waves 
with  a  slowness  vector  oiO)  the  delays  should  be  set  equal  to 

T„  =  -aoX„.  (10) 


Suppose  a  plane  wave  s{x,t)  has  a  slowness  vector  a  such  that  s(x,t)  = 
exp[j27rf{t  —  a^x)].  When  the  beam  is  steered  towards  ao,  the  output  from  the 
beamformer  due  to  5(x,t)  is  given  by  the  following: 

bfit)  =  WnTnit  -  Tn) 

n=0 
N-1 

=  TrJl  WnS{x.n,  t  +  a^X„) 

n=0 

N-1 

=  ju  I]  “'n  exp[727r/(t  +  a^x„  -  a^x„)] 

n=0 

f  I  N-1  \ 

=  I  exp[- j27r/(a  -  ao)^x„]  J  exp(j27r/t) 

=  W {f{a  -  ao))  exp(j27r/i) 

=  W(k-ko)exp(j27r/t),  (11) 


where 

=  4  Wnexp(-j27rk^x„)  (12) 

n=0 
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is  called  the  array  pattern.  The  array  pattern  indicates  how  a  signal  with  a  wavenum¬ 
ber  vector  k  =  fa  will  be  attenuated  when  the  beam  is  steered  in  the  direction  of 
ckq.  Using  the  expression  for  s(x,t)  from  Eqn.  (2)  and  the  definition  of  the  array 
pattern  from  Eqn.  (12),  the  beamformer  output  can  be  written  as 

1  /»00  /*00 

-  /“o)  exp{j2Trft)dkdf.  (13) 

iV  J—00  J—oo 

This  equation  gives  the  wideband  response  of  the  weighted  delay-and-sum  beam- 
former  and  shows  how  signals  from  different  directions  and  at  different  frequencies 
are  attenuated  before  being  combined  to  produce  the  beamformer  output. 

As  a  special  case,  consider  a  waveform  v{t),  which  may  be  a  single  signal  or 
combination  of  signals  propagating  in  the  same  direction.  Then  s(x,  t)  can  be  written 
as 

s(x,  t)  =  v{t  —  a^x)  (14) 

and  its  wavenumber-frequency  spectrum  is  given  by 

S(k,f}  =  V{f)S(k-fc,),  (15) 

where  V{f)  is  the  Fourier  transform  of  v{t)  and  (5(k)  is  the  3-D  impulse  function. 
Then  Eqn.  (13)  becomes 

/OO 

V if)W (/(a  -  ao))  exp{j27vft)df.  (16) 

-OO 

If  the  direction  of  propagation  of  v{t)  coincides  with  the  direction  to  which  the 
beamformer  is  being  steered,  then  a  =  ao  and  the  signal  waveform  will  not  be 
attenuated. 

If  q:  /  aO)  tffe  term  /(a  —  Oq)  grows  linearly  with  frequency,  causing  the 
higher-frequency  components  to  be  attenuated  more  than  lower-frequency  compo¬ 
nents.  This  effect  can  be  seen  in  Fig.  4.  Signals  from  the  ko  direction  will  be  passed 
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Passband  of  W(k  —  ko) 


Figure  4  Lowpass  Filtering  of  a  misaligned  signal. 

unattenuated  at  all  frequencies,  while  signals  originating  from  other  directions  will 
have  their  high-frequency  components  filtered  out. 

It  is  possible  to  interpret  Eqn.  (13)  as  a  4-D  filtering  operation.  Applying  the 
convolution  theorem  to  Eqn.  (7)  yields 

1  roo  /•oo  „ 

^  Loo  Loo  -  k^x)]dkd/.  (17) 

Equation  (13)  is  equivalent  to  Eqn.  (17)  evaluated  at  k  =  0.  Therefore,  the 
wavenumber-frequency  response  is  given  by 

H{kJ)  =  W{k-fao).  (18) 

Evaluating  the  array  pattern  evaluated  at  k  -  fao  will  produce  the  complex  ampli¬ 
tude  of  the  response  of  the  filter  used  to  process  the  space-time  signal  s(x,t). 

The  simplest  type  of  antenna  array  consists  of  N  identical  receiver  elements 
arranged  in  a  line  with  equal  spacing  d  between  elements.  This  configuration  is 
called  a  uniform  linear  array.  For  now,  equal  weights  Wn  =  1  will  be  assumed 
for  each  element.  If  the  array  is  oriented  along  the  x-axis,  the  position  of  the  nth 
element  is  given  by  x„  =  (nd,  0, 0)  (n  =  0, 1, . , . ,  AT  —  1)  and  the  expression  for  the 
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Figure  5  Array  pattern  for  a  10  element  linear  array  with  no  time  delays. 


array  pattern  in  Eqn.  (12)  becomes 


=  4  S  exp[-j27rk^x„] 

n=0 


N-1 


N 


N 


exp[-j27rkxnd\ 

n=0 

^  X)  (exp[-j27rA:a;d])" 
n=0 

1  1  —  exp[— j27riV’A:a;(i] 

N  1  —  exp[— j27rfea;<i] 

1  exp[j7riVfcj:d]  —  exp[—j7rNkxd]  exp[—jTTNkxd\ 
N  exp[jnkxd]  —  exp[— jf7rfca;d]  exp[—j7vkxd] 
siii{7rNkxd) 


N  sm{7rkxd) 


exp[-j7r(iV  -  l)kxd\, 


(19) 


where  kx  =  |k|  sin(0o)  =  sin(0o)/A  and  (f>Q  is  as  shown  in  Fig.  4.  The  magnitude  of 
this  array  pattern  is  shown  in  Fig.  5  for  iV  =  10.  The  pattern  repeats  with  period 
kx  =  1/d.  The  repeated  main  lobes  are  called  grating  lobes.  These  grating  lobes  can 
allow  signals  to  enter  the  beamformer  from  directions  and  frequencies  other  than 
those  for  which  the  beamformer  is  designed.  This  eflPect  is  illustrated  in  Fig.  6. 
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Figure  6  Effect  of  grating  lobes  in  (k,  /)-space.  Components  of  kx  which  lie  within 
the  grating  lobes  will  enter  the  beamformer  unattenuated. 

A  signal  with  a  wavenumber  vector  fcoi  =  fotQx  will  lie  entirely  within  the  main 
lobe.  A  signal  with  wavenumber  vector  kx  =  fot  will  cross  the  grating  lobes,  and 
components  of  this  signal  which  lie  within  the  main  lobe  or  grating  lobes  will  enter 
the  beamformer  unattenuated.  The  existence  of  grating  lobes  indicates  the  spatial 
sampling  interval  d  is  too  large  in  relation  to  the  wavelengths  of  the  incoming  signals. 
This  effect  is  equivalent  to  the  aliasing  which  results  when  time  varying  signals  are 
not  sampled  at  a  high  enough  rate.  To  avoid  aliasing  when  sampling  time  varying 
signals,  the  sample  rate  must  be  at  least  twice  the  maximum  frequency  contained  in 
the  signal.  Similarly,  to  avoid  the  spatial  frequency  aliasing  which  results  in  grating 
lobes,  the  interelement  distance  d  must  meet  the  criterion 

d  <  (20) 

where  Amin  is  the  minimum  wavelength  in  the  received  signal. 

2.3.2  Filter-and-Sum  Beamforming.  The  weighted  delay-and-sum  beam- 
former  discussed  above  can  be  extended  to  allow  the  receiver  weights  to  be  a  function 
of  frequency,  so  that  signal  components  at  different  frequencies  are  weighted  differ- 
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ently.  Then,  at  a  specific  frequency,  Eqn.  (9)  can  be  rewritten  as 


W./)  =  4EV„(/)r„((-T„),  (21) 

n=0 

where  r„(t)  is  the  inverse  Fourier  transform  of  the  signal  spectrum  incident  on  the 
nth  receiver, 

/oo 

Rn{f)exp{j2TTft)df.  (22) 

■OO 

For  a  single  frequency  /  =  /o,  Eqn.  (22)  reduces  to 


rn{t)  =  Rnifo)  exp(;27r/ot),  (23) 

and  Eqn.  (21)  can  be  rewritten  as 

bf{t,  fo)  =  ^Yl  Wn{fQ)Rnifo)  exp[j2nfo{t  -  Tn)].  (24) 

Integrating  this  expression  over  all  frequencies  gives  the  filter-and-sum  beamformer 
output 

Mt)  =  rbf{t,f)df 

00 

=  Wn{f)Rn{f)exp[j2nf{t-rn)]f.  (25) 

To  simplify  the  notation,  define  the  signal 

qn{t)=  f  W„{f)Rn{f)exi>[j2'Kft]df.  (26) 

J—00 

Then  Eqn.  (25)  can  be  written  as 

f^i*)  =  (27) 

n=0 
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Prom  Eqn.  (26),  it  can  be  seen  that  g„(t)  is  the  convolution  of  Wn{t)  with 
where  Wn{t)  is  the  inverse  Fourier  transform  of  the  weighing  function  Wn{f).  The 
function  Wn{t)  can  be  viewed  as  the  impulse  response  of  a  filter  that  operates  on  the 
receiver  signal  The  filtered  signal  qn{t)  which  results  is  then  used  to  shape  the 
beam. 


2.3.3  Frequency- Domain  Beamforming.  It  is  possible  to  write  expressions 
for  the  outputs  of  the  weighted  delay-and-sum  and  filter-and-sum  beamformers  in 
the  frequency  domain  by  computing  their  Fourier  transform  as  follows 

BF{f)  =  j-'^WnRn{f)exp(-j2TrfTn)  (28) 

n=0 

(29) 

•'''  n=0 

The  output  of  each  beamformer  at  a  specific  frequency  /  is  obtained  by  multiply¬ 
ing  the  above  expressions  by  exp{j2iTft),  and  the  complete  beamformer  output  is 
obtained  by  integrating  over  all  frequencies 


II 

f  BF{f)  exp{j2nft)df 
«/-00 

(30) 

fs{t)  = 

roo 

1  FS{f)exp{j2Trft)df. 

J—OO 

(31) 

Since  computing  the  signal  spectrum  Rn{f)  in  Eqns.  (28)  and  (29)  would  re¬ 
quire  integrating  the  signal  over  all  time,  it  is  necessary  to  approximate  the  spectrum 
by  computing  the  Fourier  transform  of  a  segment  of  the  time  signal.  This  is  done 
by  applying  a  finite-extent  window  to  the  signal 

PT{t-T)rn{r)exp{-j2TTfT)dT,  (32) 

J—OO 

where  pxit  —  r)  is  the  window  function  centered  at  time  r  with  width  T. 
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Using  this  definition  in  Eqn.  (28)  gives  the  frequency-domain  beamformer  out¬ 
put  of  the  weighted  delay-and-sum  beamformer  as 

/)  =  4  f)  exp[j27r/(t  -  r„)].  (33) 

n=0 

Integrating  this  expression  over  all  frequencies  will  give  an  approximation  to  the 
beamformer  output  signal.  A  similar  expression  results  if  Eqn.  (29)  is  used. 

Using  Eqns.  (1),  (10),  (12)  and  defining 

r„(<)  =  s(x,i,t)  (34) 

P(f)  =  /  PTit)exp{-j2'Kft)dt,  (35) 

J—OO 

it  is  possible  to  write  the  frequency  domain  beamformer  output  as 

1  roo  /•oo 

/-»  (36) 

If  the  array  pattern  lU(k  —  /ao)  is  considered  to  be  zero  except  near  k  =  fcxo  and 
the  window  spectrum  P{9  —  f)  is  considered  to  be  zero  except  near  B  =  f,  the  factor 
W(k  —  fao)P{B  —  f)  will  represent  the  passband  of  the  system.  The  beamformer 
will  attenuate  all  signals  in  (k,  0)-space  except  those  near  k  =  /ao  and  9  =  f.  The 
passband  is  illustrated  in  Fig.  7. 

2.4  Discrete  Time  Beamforming 

Until  this  point,  all  signals  have  been  treated  as  continuous  time  signals.  How¬ 
ever,  most  beamformers  are  implemented  digitally,  requiring  the  received  signals  to 
be  sampled  in  time.  In  this  section,  the  effects  of  using  discrete  time  signals  will  be 
examined.  The  same  notation  will  be  used  as  with  the  continuous  time  case,  but 
now  the  variable  t  will  represent  an  integer  multiple  k  of  the  sampling  period  T,  such 
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Figure  7  Bandpass  filter  in  wavenumber-firequency  space.  The  shaded  region  rep¬ 
resents  the  passband. 


that  t  =  kT.  Thus,  r„(t)  will  denote  a  discrete-time  signal.  It  is  assumed  that  the 
sampling  rate  is  high  enough  to  avoid  aliasing. 

The  discrete-time  version  of  the  weighted  delay-and-sum  beamformer  is  given 
by 

l>f{kT)  =  ^'£w„r„{kT-k„T),  (37) 

n=0 

where  knT  is  the  quantized  steering  delay  for  the  nth  sensor.  Since  is  an  integer, 
the  quantized  steering  delays  will  generally  not  be  equal  to  the  ideal  delay  Tn-  The 
quantization  error  in  the  nth  delay  is  given  by 


At,i  —  kjiT  Tfi- 


This  error  must  be  accounted  for  in  the  response  of  the  beamformer  to  a  plane  wave 
with  slowness  vector  ao.  The  new  expression  is 

^  N-l 

f)  =  Tr'^'^n  exp[-j27r(k  -  /ao)^x„]  exp(-j27r/<Arn).  (39) 

n=0 
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Figure  8  Quantized  Steering  Delays.  Dots  indicate  the  quantized  steering  delays, 
the  line  indicates  the  ideal  delay. 


As  a  simple  example  to  illustrate  the  effect  of  quantizing  error,  consider  a  uniform 
linear  array  oriented  along  the  rc-axis.  If  the  beam  is  steered  such  that  ao®  =  T/2d, 
then  the  ideal  delay  is  =  —nT/2.  However,  the  delay  applied  must  be  an  integer 
multiple  of  the  sample  period.  Therefore,  the  actual  delay  will  be  knT,  where  is 
the  largest  integer  less  than  or  equal  to  n/2.  The  quantized  steering  error  will  be 


=  < 


^  for  n  odd 
0  for  n  even 


(40) 


as  shown  in  Fig.  8.  For  this  case,  the  frequency  response  can  be  written  as 


,  K/2-1 
^  ?o 

I  N/2-1 

+  ]v  E  exp 


n=0 


-j2'!r  f fcx  -  {2n  +  l)d  - 


f  1  1 

j  E  exp[-j27rn(2A:a:d  -  ft)]  HI  +  exp{-j27rk^d)] 


=  H,{kJ)H2{k,f), 


(41) 
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Figure  9  Effect  of  steering  delay  quantization  on  array  pattern:  N  =  8,  d  =  1,T  = 
1,  /  =  1/4.  (a)  Quantized  beamformer  [Eqn.  (41)],  (b)  Ideal  beamformer 
[Eqn.  (19)]. 

where 

An  example  of  the  magnitude  (k,  /)  |  is  shown  in  Fig.  9  for  both  the  quantized 
and  ideal  cases.  For  the  quantized  case,  the  main  lobe  peak  is  lower  and  some  of 
the  side  lobes  are  higher  than  for  the  ideal  case.  The  ideal  response  is  given  by 
/)!  =  \W{kxd  —  /r/2)|,  where  PF’(k)  is  as  defined  in  Eqn.  (19). 

2.5  Total  Antenna  Pattern 

To  compute  the  total  antenna  pattern  for  a  given  array,  the  array  pattern 
computed  using  Eqn.  (19)  must  be  multiplied  by  the  element  pattern.  For  a  small 
dipole,  the  element  pattern  will  be  a  cosine  function,  with  the  angle  measured  from 


2  N/2-1 

exp[-j27rn(2M-/0] 

n=0 

sin[(7riV)(M-/r/2)]  f 


(7riV)sin(M-/r/2) 
=  ^[H-exp(-j27rM)] 

=  exp{—jiTkxd)  cos{‘Kkxd). 


(f  -  l)  (, 
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the  normal  to  the  element.  Substituting  kx  =  sin(^)/A  into  the  equation  for  the 
array  pattern  in  Eqn.  (19)  and  multiplying  by  the  element  pattern,  the  magnitude 
of  total  antenna  voltage  pattern  can  be  computed  as  a  function  of  angle  as 


sin(7riVdsin((^)/A)  _ 
APJ(p)  =  .  /  cos(0). 

^  ’  N  sm(7rd  sm((^)  /A) 


(44) 


The  antenna  power  pattern  is  found  by  squaring  Eqn.  (44), 


APM 


sin(7riVd  sin(^)  /  A) 


[N  sm{ndsm{(f>)  /  X) 


cos(0) 


(46) 


The  gain  of  the  array  is  equal  to  the  inverse  of  the  average  radiation  intensity.  The 
radiation  intensity  of  the  array  is  given  by  [2:276-278] 


sin(7rJVd  sin(^)/A) 


l2 


N  sin(7rd  sin((f>)  /  A)  J 


(46) 


The  average  radiation  intensity  of  the  array,  l7o,  is  given  by 


Uo 


/•ir/2 

sin(7riVd  sin((/>)/A) 

/— ir/2 

N  sin(7rd  sin(0)/A) 

cos{<l>)d<l>. 


(47) 


By  letting  ■0  =  27rdsin(0)/A,  the  average  radiation  intensity  can  be  expressed  as 


sin(iV’V’/2) 

iV’sin('0/2) 


di’, 


(48) 


and  the  array  gain  =  I/Uq.  If  the  spacing  between  array  elements  is  equal  to 

a  half  wavelength,  the  array  gain  is  approximately  equal  to  the  number  of  elements. 
To  find  the  total  antenna  transmit  gain,  Gt{4>),  the  array  gain  is  multiplied  by  the 
element  gain,  ^(0). 
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2.6  Summary 

This  chapter  has  presented  some  of  the  fundamentals  of  beamforming  with  a 
linear  array.  Expressions  were  derived  for  the  weighted  delay-and-sum  beamformer 
and  the  filter-and-sum  beamformer  in  both  the  time  and  frequency  domains.  The 
effects  of  quantization  error  introduced  by  a  digital  implementation  of  a  weighted 
delay-and-sum  beamformer  were  also  examined.  Finally,  an  expression  for  the  gain  of 
a  linear  array  was  presented.  The  next  chapter  presents  a  theoretical  development  of 
space-time  adaptive  processing,  which  is  an  extension  of  the  beamforming  techniques 
presented  in  this  chapter.  It  will  be  assumed  there  is  no  quantization  error. 
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III.  Space- Time  Adaptive  Processing 

3.1  Introduction 

Now  that  we  have  the  basics  of  multidimensional  signal  processing,  we  will 
begin  studying  the  signals  in  more  detail  to  understand  how  they  are  composed 
and  how  they  can  be  analyzed.  In  this  chapter,  models  are  developed  for  the  radar 
system  and  signals  from  the  target,  jammers,  and  clutter.  The  chapter  concludes 
with  a  discussion  of  fully  adaptive  STAR  This  discussion  closely  follows  that  in  [8]. 

3.2  Radar  System 

The  system  being  analyzed  is  assumed  to  be  mounted  on  an  airborne  platform 
and  utilizing  a  pulsed  Doppler  waveform.  The  platform  is  at  a  height  ha  and  traveling 
with  a  constant  velocity,  Va.  The  variables  (j)  and  6  are  the  azimuth  and  elevation 
angles  measured  with  respect  to  the  normal  to  the  array.  The  unit  vector  k  indicates 
the  direction  to  the  point  of  interest  and  may  be  represented  in  Cartesian  coordinates 
as 

k(0, 0)  =  cos{9)  sin((^)x  +  cos(0)  cos(0)y  +  sin(0)z.  (49) 

Figure  10  depicts  the  scenario.  The  antenna  is  assumed  to  be  a  uniform  linear  array 
of  N  identical  elements  separated  by  a  distance  d  and  oriented  along  the  ar-axis. 
The  vector  d  =  dx  will  be  used  to  represent  the  array  spacing  and  orientation. 
Each  element  has  the  same  radiation  voltage  pattern  f{<l>,9)  and  power  pattern 
g{^,9)  =  |/(^,  0)p.  The  radar  transmits  a  coherent  burst  of  M  pulses  on  a  carrier 
frequency  /o  =  c/Aq  and  at  a  constant  pulse  repetition  frequency  (PRF),  =  l/T^, 
where  %  is  the  pulse  repetition  interval  (PRI).  The  time  interval  over  which  the 
waveform  returns  are  collected  is  equal  to  MT^  and  is  called  the  coherent  processing 
interval  (CPI). 
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Figure  10  Radar  system  and  antenna  array  geometries 

The  radar  array  has  a  transmit  pattern  Gti(f>,d)  which  is  determined  by  the 
number  and  spacing  of  the  elements.  Each  element  is  connected  to  its  own  receiver, 
consisting  of  a  down-converter,  matched  filter,  and  A/D  converter.  The  receiver 
bandwidth  B  is  taken  to  be  equal  to  the  transmitted  pulse  bandwidth.  The  output 
from  each  element’s  A/D  converter  is  sent  to  the  processor,  where  the  total  signal 
is  analyzed.  For  each  PRI,  L  time  samples  are  collected  to  cover  the  range  interval 
of  interest.  If  the  array  has  N  receivers  and  transmits  M  pulses,  the  received  data 
during  a  single  CPI  will  consist  of  LMN  complex  samples.  The  data  can  be  conve¬ 
niently  represented  as  &  L  x  M  x  N  cube  of  complex  samples  referred  to  as  the  CPI 
datacube,  as  illustrated  in  Fig.  11. 

For  the  Zth  range  gate,  Xm,!  is  an  iV  x  1  vector  of  all  the  receiver  element 
outputs  for  the  mth  pulse.  This  vector  is  called  a  spatial  snapshot.  The  matrix 

Xj  is  an  A'  X  M  matrix  consisting  of  the  spatial  snapshots  from  all  pulses.  By 

• 

stacking  the  columns  of  Xj  on  top  of  each  other,  this  matrix  is  reshaped  into  an 
MN  X  1  vector  Xi  =  vec(Xi)  called  the  space-time  snapshot  The  remainder  of  the 
analysis  will  assume  the  returns  from  only  one  range  gate  are  being  studied,  and 
the  I  subscript  will  be  dropped.  Therefore,  x  will  be  the  space-time  snapshot  at  the 
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Figure  11  Radar  CPI  datacube 

range  of  interest,  and  will  represent  the  spatial  snapshot  from  the  mth  PRI  at 
this  range. 

The  space-time  snapshot  has  two  components,  the  target  return  and  the  un¬ 
desired  signal, 

X  =  Xt  +  Xu  =  Oit^t  +  Xu)  (50) 

where  at  is  the  complex  amplitude  of  the  target  return,  Vf  is  the  system  response  to 
a  target  of  unit  amplitude,  and  Xu  =  Xc  +  Xj  +  Xn  represents  undesired  signals,  due 
to  clutter  return,  jamming  and  noise.  If  no  target  is  present,  at  =  0  and  the  space- 
time  snapshot  consists  entirely  of  the  undesired  signals.  The  presence  of  a  target 
causes  a  shift  in  the  mean  of  the  data,  but  does  not  affect  the  covariance  of  the  data. 
Assuming  the  clutter,  jamming  and  noise  signals  to  be  mutually  uncorrelated,  the 
covariance  matrix  of  the  total  undesired  signal  is  the  sum  of  the  covariance  matrices 
of  each  component: 

Ru  =  -E^[XuX^]  =  ^[(Xc  +  Xj  +  Xj(Xc  +  Xj  +  Xn)^] 

=  E[XcXc]  +  E[XjXf]  +  E[XnXn] 

=  Rc  +  Rj  +  Rn)  (51) 
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where  the  subscripts  c,  j,  and  n  indicate  clutter,  jammer,  and  noise  interference 
respectively.  The  following  sections  develop  a  model  for  each  component  of  the 
space- time  snapshot. 


3.3  Target 

The  return  from  a  target  at  a  given  range,  Rt,  is  a  function  of  the  target’s 
azimuth  angle,  0,  depression  angle,  6,  closing  velocity  with  respect  to  the  radar,  Vt, 
and  radar  cross-section  (RCS),  df  The  received  signal  can  be  derived  by  defining 
the  transmitted  signal  and  applying  the  filtering  that  produces  the  samples  of  the 
space- time  snapshot.  The  transmitted  coherent  burst  of  pulses  from  the  array  is 
represented  by 

s(t)  =  ,  (52) 


where 


M-l 


m=0 


(53) 


is  the  complex  envelope  of  the  signal,  is  the  PRI,  and  Up{t)  is  the  complex  envelope 
of  a  single  pulse  of  width  Tp  and  having  unit  energy.  The  transmitted  signal  has 
amplitude  at  and  a  random  phase  ijj.  The  signal  received  from  the  target  by  the  nth 
element  is  given  by 


Sn{t)  =  aru{t  -  r„)  exp[i27r(/o  -h  ft){t  -  r„)]e^’^,  (54) 


where  is  the  amplitude  of  the  returned  signal  and  ft  =  2vt/XQ  is  the  target  Doppler 
frequency.  The  target  delay  to  the  nth  element  has  two  components:  r„  =  Tj  -I-  r^, 
where  r*  =  2Rt/c  is  the  time  for  the  signal  to  travel  from  the  array  to  the  target  and 
back,  and 


,  nk(0t,0t)-d„  d  ^  , 

r„  = - =  — n-  cos  0t  sin  (pt 

c  c 


(55) 
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Figure  12  Illustration  of  relationship  between  wavenumber  vector  k,  spatial  fre¬ 
quency  1?,  wavelength  A,  and  element  spacing  d. 

is  the  relative  delay  measured  from  the  phase  reference  to  the  nth  element.  The 
target  spatial  frequency  will  be  defined  as 

=  A  cos  9t  sin  (56) 

Ao  Ao 

Figure  12  shows  how  the  spatial  frequency  relates  to  the  wavelength,  array  spacing, 
and  wavenumber  vector.  The  phase  delay  to  the  nth  element  can  now  be  given  by 
-Uqt^  =  n2'K'dt.  Assuming  a  narrowband  transmitted  waveform,  the  relative  delay 
term  t!^  in  the  envelope  u{t  -  Tn)  of  Eqn.  (54)  is  insignificant.  The  received  signal 
can  be  expressed  as 

Sn(t)  =  aru(t  -  Tt)e^''^  expU27r(fot  +  fit  -  foTn  -  ftT„)] 

=  aru(t  -  Tt)e^'^  exp[j27r(/o^  +  fti  -  hn  -  foTn  -  fm  -  /t?-^)] 

=  (57) 

where  V*  now  contains  several  of  the  fixed  frequency  terms.  The  signal  is  downcon- 
verted  to  baseband,  removing  the  carrier  frequency  term,  and  then  passed  through 
a  matched  filter  with  a  response  given  by  h{t)  =  u*{-t).  The  output  of  the  matched 
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filter  for  the  nth  channel  is  found  from  the  convolution  of  the  signal  with  the  filter, 

/oo 

aru{(3  -  -  t)dl3 

-OO 

roo  A/-1 

=  ^  np()0  -  n  -  mr,)n;(/3  -  (58) 

m=0 

Letting  7  =  fi—Tt—mTr  and  interchanging  the  integration  and  summation  operations 
gives 

M-l  foo 

Xn(t)  =  5;  /  np(7)np(7-<  +  Tt  +  mT,)e^'2’^-^‘('^+"‘+”^^^)d7 

m=0 

A^-1  /-oo 

=  ^  /  np(7)n!(7  - 1  +  r*  +  mTr)e^'’"’^'^‘T'd7. 

m=0 

(59) 

The  term  represents  a  fixed  frequency  and  can  be  absorbed  into  the  random 

phase  term.  The  integral  can  be  recognized  as  the  ambiguity  function  [7:411-412], 

X{r,  f)=  f  «p(7)Mp(7  -  r)e^^^f'^d'r.  (60) 

J—00 

At  this  time  it  will  be  convenient  to  normalize  the  target  Doppler  frequency  to  the 
PRF.  The  normalized  Doppler  frequency  is  defined  as 

Jr 

With  these  two  definitions,  Eqn.  (59)  can  be  written  as 

M-l 

Xnit)  =  ^  -  r*  -  mT,,  ft).  (62) 

m—O 
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For  the  range  gate  containing  the  target,  tm  =  Tt  +  mTj.{m  =  0, —  1) 
represents  the  sample  times  from  each  PRJ  at  this  range  gate.  The  target  samples 
are  given  by 

/,).  (63) 

Because  the  pulse  waveform  is  normalized,  x(0, 0)  =  1.  It  is  also  assumed  the  pulse 
waveform  is  insensitive  to  the  target  Doppler,  such  that  x(0, /)  «  1.  Therefore, 
Eqn,  (63)  can  be  written  as 


^mn 


gim27rrot  ^ 


n=0, 1, . . . ,  AT  —  1 

m=0, 1 


(64) 


where  at  =  is  the  complex  random  amplitude  of  the  signal.  The  target  ampli¬ 
tude  can  be  calculated  from  the  radar  range  equation  as  follows.  The  single  pulse 
signal-to-noise  ratio  (SNR)  for  a  single  element  is  given  by 


(65) 


where  the  terms  are  as  defined  in  Table  1.  The  target  power  is  computed  from 


=  (66) 

where  cr^  is  the  thermal  noise  power  per  element.  The  target  amplitude  is  then  given 
by  flr  =  v^- 

At  this  time,  it  will  be  convenient  to  define  the  spatial  and  temporal  steering 
vectors.  The  spatial  steering  vector  is  an  AT  x  1  vector  given  by 


a(^,«) 


1 .  fJ2ir4-  cos 9 sin 


•  cos0sin0j 


(67) 
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Table  1  Symbols  for  Radar  System  Parameters 


N 

number  of  elements  in  array 

M 

number  of  pulses  in  CPI 

L 

number  of  range  samples  per  PRI 

fo 

radar  operating  frequency 

Ao 

radar  operating  wavelength 

d 

interelement  spacing 

Tr 

pulse  repetition  interval  (PRI) 

fr 

pulse  repetition  frequency  (PRF) 

Pt 

peak  transmit  power 

Tp 

transmit  pulse  width 

B 

instantaneous  bandwidth 

GtiOA) 

full  array  transmit  power  gain 

m<i>) 

element  pattern  (voltage) 

9{d,  <f>) 

element  pattern  (power) 

Ijf 

system  losses  on  receive 

Lt 

system  losses  on  transmit 

Ls 

total  system  losses 

iVo 

receiver  noise  power  spectral  density 

Rt 

target  range 

target  RCS 

h(i 

platform  height 

Va 

platform  velocity 
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or 


(68) 


The  temporal  steering  vector  b(t<7)  is  an  Af  x  1  vector  given  by 

b(n7)  =  [l;  . . . ;  ,  (69) 

where  d  and  vj  are  as  defined  in  Eqns.  (56)  and  (61).  The  MN  x  1  space-time 
steering  vector  will  be  defined  as 

V  =  v(i?,  w)  =  b(ti7)  (g)  a(i9),  (70) 

where  ((g)  indicates  the  Kronecker  product.  With  these  conventions,  the  space-time 
snapshot  of  the  target  data  can  be  written  as 


Xt  =  Oith{wt)  (8>  a(i?f)  =  Q!tb<  0  a<  =  Q;<Vt,  (71) 

which  is  the  target  component  of  the  total  space- time  snapshot  of  Eqn.  (50).  The 
space-time  steering  vector  of  the  target  Vf  is  also  referred  to  as  the  target  steering 
vector. 


3.4  Noise 

Noise  is  be  present  in  each  of  the  receiver  elements.  Since  each  element  has 
its  own  receiver,  it  is  assumed  the  noise  processes  on  each  of  the  N  element  will  be 
mutually  uncorrelated.  It  is  also  be  assumed  that  the  M  noise  samples  on  a  given 
element  separated  by  a  nonzero  multiple  of  the  PRI  are  temporally  uncorrelated. 
These  two  assumptions  may  be  expressed  as 


^  {®ni,m2?n2,Tn}  ~  <^n^ni-n2 
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(72) 

(73) 


where  =  NqB  is  the  noise  power  per  element  and  Sm  is  the  Kronecker  delta 
function  defined  as 

1  ,  m  =  0  ,  ^ 

(74) 

0  ,  m^O 

The  noise  component  of  the  space-time  covariance  matrix  can  then  be  expressed  as 

=  E  {XnXn  }  =  0  IjV  =  (75) 

The  noise  power  will  be  normalized  to  =  1  so  that  all  signal  levels  can  be  refer¬ 
enced  by  their  SNR  per  element  and  pulse. 

3.5  Jamming 

For  this  analysis,  only  barrage  noise  jamming  from  a  source  at  a  long  range  and 
an  angle  different  from  that  of  the  target  will  be  considered.  The  jammer  bandwidth 
will  be  assumed  to  be  much  larger  than  the  pulse  bandwidth,  and  the  radar  PRF 
will  be  assumed  to  be  much  smaller  than  the  instantaneous  bandwidth.  The  jammer 
signal  will  be  correlated  spatially  across  the  array,  but  temporally  uncorrelated  from 
pulse  to  pulse.  Therefore,  the  jammer  signal  will  appear  like  thermal  noise  in  the 
time  domain,  but  as  a  point  source  in  the  spatial  domain. 

For  a  jammer  at  range  Rj,  elevation  $j,  and  azimuth  <pj  with  effective  radiated 
power  density  Sj,  the  jammer  power  spectral  density  Jq  received  by  one  element  is 
given  by 

_  Sjg{9j,  <^j)Ao 

{4nyRjLr  ' 

The  jammer-to-noise  ratio  (JNR)  at  the  element  is  given  by  =  Jq/Nq.  The  jammer 
space-time  snapshot  is  given  by 


Xj  =  ctj  ®  a^, 


(77) 
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where  aj  is  a  M  x  1  vector  of  the  random  jammer  amplitudes  during  each  PRI 
and  aj  =  $j)  is  the  jammer  steering  vector.  Assuming  the  jamming  signal  is 

stationary  over  a  CPI,  the  jammer  space-time  covariance  matrix  can  be  calculated 
as 


^3=^{XjXf}  =  ^?{(aj0aj)(Q:j0aj)^} 

=  £?{(Q!j0aj)(af  0af)} 
=  (aj-af)} 

=  <S)  {ajaf) 

=  1m  <B>  (crlCj^jaf) 

=  1m® 


(78) 


where  is  the  jammer  spatial  covariance  matrix. 

This  result  can  be  generalized  to  multiple  jamming  sources.  If  there  are  J 
jammer  sources,  the  spatial  covariance  matrix  is  given  by 

=  AjSjAf,  (79) 

where 

Aj  =  [a(0i0i),  a(0202),-.-,  a(0j0j)]  (80) 

is  the  N  X  J  jammer  spatial  steering  matrix.  If  the  jammer  waveforms  are  mutually 
uncorrelated,  the  jammer  source  covariance  matrix  Sj  will  be  a  J  x  J  diagonal  matrix 
of  the  jammer  powers.  Prom  Eqn.  (78)  it  can  be  seen  that  the  jammer  space-time 
covariance  matrix  will  be  block  diagonal,  with  each  non-zero  block  consisting  of  the 
jammer  spatial  covariance  matrix.  If  the  J  jammer  steering  vectors  are  linearly 
independent  and  no  two  jamming  signals  are  perfectly  coherent,  then  the  rank  of 
will  be  equal  to  J,  and  Rj  will  have  rank  MJ. 
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Radar  Platform 


Figure  13  Clutter  Geometry 


3.6  Clutter 

Clutter  is  defined  as  the  radar  return  from  the  ground.  It  is  the  most  compli¬ 
cated  source  of  interference  because  it  is  distributed  in  angle,  range,  and  Doppler. 
It  exists  over  all  azimuth  angles  and  all  elevation  angles  below  the  horizon.  The 
clutter  range  returns  extend  from  the  platform  altitude  to  the  radar  horizon,  and 
the  Doppler  returns  cover  a  range  of  frequencies  corresponding  to  —Va  to  -\-Va. 

Figure  13  depicts  a  clutter  ring  at  a  range  Rc.  A  4/3  effective  radius  is  assumed 
to  model  the  curvature  of  the  earth,  so  the  effective  earth  radius  Og  =  (4/3)re.  The 
elevation  angle  to  the  clutter  ring  is  given  by 


Q  a  /T>  \  ’-if  ^a{ha  +  20e)^  --if  ,  ha 

0.  =  W)  =  -8m  (  2R.{a.  +  K)  +  h:, 


(81) 


and  the  grazing  angle  is  given  by 


I  I  f  TJ  \  '—if  ha(ha  +  20e)^  .  f  Rc 

=  =  ( - — - l«-8m  (--- 


(82) 


The  grazing  angle  is  the  angle  between  the  vector  from  the  antenna  to  the  clutter 
patch  and  the  tangent  at  the  clutter  patch. 
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The  radar  horizon  range  is  given  by  Rh  =  y^2aeho  +  hi  and  the  radar’s  unam¬ 
biguous  range  is  given  by  =  c/2fr.  li  Ru  >  Rh,  the  clutter  will  be  unambiguous 
in  range  and  the  return  in  the  Ith.  range  gate  will  contain  clutter  from  at  most  one 
range.  If  Ru  <  Rh,  the  clutter  will  be  range  ambiguous  and  some  or  all  of  the  range 
gates  will  contain  clutter  returns  from  multiple  ranges.  At  this  time,  the  assumption 
will  be  made  that  the  clutter  is  unambiguous  in  range,  and  that  the  clutter  return 
in  a  particular  range  gate  is  from  a  single  clutter  ring.  The  clutter  return  will  be 
modeled  as  the  superposition  of  Nc  independent  clutter  sources  evenly  distributed 
in  azimuth. 

The  location  of  a  particular  clutter  patch  will  be  described  by  its  azimuth  angle 
and  its  range  Rc  (or  elevation  angle  9c).  The  spatial  frequency  for  the  kth.  clutter 
patch  is  given  by 

.  d  „  .  ^ 

Vk  =  — — —  =  -p  COS  9c  sia  4>k  (83) 

Aq  Aq 

and  the  normalized  Doppler  frequency  is  represented  by  ot*.  The  clutter  component 
of  the  space-time  snapshot  is  given  by 

Nc  Nc 

Xc  =  £  2  <S»  Bfc,  (84) 

fc=i  *=1 

where  is  the  complex  random  amplitude  from  the  fcth  clutter  patch  and  vj,  = 
v(^?ft,it7fe)  is  the  space-time  steering  vector  to  the  fcth  clutter  patch,  as  defined  in 
Eqn.  (70). 

Each  clutter  patch  represents  an  effective  area  bounded  by  the  radar’s  reso¬ 
lution  in  azimuth,  =  2Tr /Nc,  and  range,  AR  =  c/2B.  The  effective  RCS  of  a 
clutter  patch  at  range  Rc  is  given  by 

<^k  =  o’o(^c)  ^k)  X  PatchArea  =  ao{9c,  (pk)Rc  AR  sec  tl^c,  (85) 
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where  ao  is  the  area  reflectivity  of  the  ground  at  the  A:th  clutter  patch.  It  will  be 
assumed  that  the  reflectivity  at  a  given  range  is  uniform  in  azimuth,  so  uo  =  7  sin 
where  7  depends  on  the  terrain.  With  these  deflnitions,  the  clutter-to-noise  ratio 
(CNR)  is  defined  as 

p  ^  PtTpGt{6c,  (f>k)9{0c, 

{4ir)moLsRt 

The  amplitude  of  the  clutter  from  the  feth  clutter  patch  is  given  by  £?  {|ajfcp}  = 

If  the  returns  from  different  clutter  patches  are  uncorrelated,  then 


£:Kar)=<^Sef5»-i-  (87) 

Using  Eqns.  (84)  and  (87),  the  space- time  covariance  matrix  for  the  clutter  is  given 
by 


=  0  a^af ,  (88) 

*=1 

where  b*  =  b(t37fc)  and  a*  =  a(i?fc).  The  clutter  spatial  and  temporal  covariance 
matrices  are  coupled  because  the  clutter  Doppler  is  a  function  of  angle.  The  clutter 
space-time  covariance  matrix  can  also  be  expressed  as  a  product  of  matrices: 

R,  =  VeBeV»,  (89) 

where  Vc  =  [vi,  V2, . . . ,  v^y,,]  is  an  MN  x  Nc  matrix  of  the  clutter  space-time  steering 
vectors,  and 

Sc  =  a^diag[Ci,  I2,  ■  •  • ,  Cn,]  (90) 

contains  the  clutter  power  distributions.  Due  to  the  assumptions  of  an  ideal  uniform 
linear  array  and  mutually  uncorrelated  clutter  patches,  the  clutter  space-time  co- 
variance  matrix  has  a  special  structure  referred  to  as  Toeplitz-block-Toeplitz.  This 
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structure  can  be  seen  in  Fig.  14.  All  blocks  along  a  diagonal  are  identical,  and  within 
each  block,  all  elements  along  a  diagonal  are  equal.  The  blocks  correspond  to  the  M 
pulses,  and  the  elements  within  each  block  correspond  to  theiV  antenna  elements. 

3.6.1  Effects  of  Platform  Motion  on  Clutter.  The  above  derivation  as¬ 
sumed  the  array  was  stationary  while  the  data  was  being  collected.  When  the 
array  is  moving,  as  it  would  be  on  an  airborne  platform,  the  clutter  develops  a 
special  structure  due  to  the  dependence  of  the  Doppler  frequency  on  angle.  As¬ 
sume  the  platform  is  moving  such  that  the  velocity  vector  of  the  array  is  given  by 
Va  =  ■ya(cos  +  sin  <l>ay),  where  (j)a  is  the  misalignment  angle  between  the  velocity 
vector  of  the  array  and  the  array  axis,  as  shown  in  Fig.  10.  This  may  occur  when  the 
aircraft  is  “crabbing”  to  compensate  for  a  crosswind.  The  clutter  Doppler  frequency 
of  a  specific  clutter  patch  is  then  given  by 

041 

h  +  (l>a)  =  -^  COS  Oc  sin(0fc  -f-  <^o)-  (91) 

The  clutter  spatial  frequency  and  normalized  Doppler  will  be,  respectively, 

'dk  =  vsm<f>k  (92) 

zz7jfc  =  I3rism{(pk  +  <l>a)  (93) 

where  =  2VaTj./d  and  t]  =  dcosOc/Xo. 

Equations  (92)  and  (93)  can  be  shown  to  satisfy  the  standard  equation  of  an 
ellipse  in  (i?,  zj7)-space, 

^2^2  _  ^2^  cos  0o)^?c^c  +  ~  <f>a  =  0.  (94) 

The  elliptical  shape  of  the  clutter  results  from  returns  entering  the  receiver  through 
the  back  lobe  of  the  antenna.  The  ellipse  will  be  rotated  clockwise  from  the  positive 
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Figure  14  Clutter  covariance  matrix  (^  =  1),  showing  the  Toeplitz-block-Toeplitz 
structure. 
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zz7-axis  by  an  angle  such  that  [8] 


cot  2^a  = 


1-13^ 
2(3  cos  ((>a 


(95) 


and  will  have  semimajor  and  semiminor  axes  given  respectively  by 


a  =  sin <f>a  (^^{(3^  +  1)  +  |c|) 
b  =  pr)sin(l)a  +  1)  -  |c|) 


-1/2 

-1/2 


where  c  is  given  by 


(96) 

(97) 


c  =  —  l)cos2^a  —  ^cos0osm2^tt. 


(98) 


In  the  case  where  the  velocity  vector  of  the  platform  is  aligned  with  the  axis 
of  the  array,  4>a  =  Oj  and  Eqn.  (94)  reduces  to  the  equation  of  a  line  of  clutter. 


'^e  =  (99) 

The  back  lobe  return  is  still  present,  but  directly  overlays  the  front  lobe  return. 
When  (3  <  1,  there  will  be  a  range  of  normalized  Doppler  frequencies  which  will  be 
free  from  clutter  return.  For  values  of  >  1,  the  clutter  return  will  be  ambiguous 
in  Doppler.  Figures  15  through  17  show  the  plots  of  the  clutter  ridges  for  various 
values  of  (3  and  different  misalignment  angles  (j>a  as  functions  of  the  spatial  frequency 
and  normalized  Doppler  frequency. 

3.6.2  Rank  of  the  Clutter  Covariance  Matrix.  The  severity  of  the  clutter 
scenario  and  the  number  of  degrees  of  freedom  needed  to  effectively  cancel  the  clutter 
are  indicated  by  the  rank  of  the  clutter  covariance  matrix.  The  rank,  rc,  can  be 
estimated  from  Brennan’s  rule: 
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Spatial  Frequency  Spatial  Frequency 

(c)  (d) 


Figure  15  Clutter  ridges  for  different  values  of  /3,  (j>a  =  0,  (a)  /5  =  0.5,  (b)  13  =  1,0, 
(c)  l3  =  1.5,  (d)  j3  =  2.8.  (The  solid  line  represents  the  front  lobe  and 
the  dashed  line  represents  the  back  lobe.) 
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Figure  16  Clutter  ridges  for  different  misalignment  angles,  /?  =  1,  (a)  =  0°, 

(b)  =  10°)  (c)  0o  =  45°,  (d)  (j)a  =  90°. (The  solid  line  represents  the 

front  lobe  and  the  dashed  line  represents  the  back  lobe.) 


Figure  17  Clutter  ridges  for  different  combinations  of  (j)a  and  /?. 
(a)  p  =  2.8,  <l>a  =  10°,  (b)  P  =  2.8,  <j)a  =  60° 
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Figure  18  Clutter  eigenspectra.  Values  of  /?  are  indicated  at  the  top. 

Theorem  1  (Brennan’s  Rule)  The  rank  of  the  clutter  covariance  matrix  is  ap¬ 
proximately  given  by 

rc  «  +  (M  -  1)/?J  (100) 

where  f3  =  2vaTr/d  is  the  number  of  half -interelement- spacings  traveled  by  the  plat¬ 
form  during  one  PRL  (The  brackets  ([J)  indicate  rounding  down  to  the  nearest 
integer.) 

When  the  eigenvalues  of  the  clutter  covariance  matrix  are  calculated,  rc  predicts 
where  the  plot  of  the  eigenspectrum  will  begin  to  roll  off.  The  sharpness  of  the 
roll-off  depends  on  the  platform  speed  Vg,,  misalignment  angle  0a,  back  lobe  level, 
and  the  existence  of  intrinsic  clutter  motion  (see  Section  3.6.3).  A  plot  of  the  eigen¬ 
spectrum  for  various  values  of  /?  is  shown  in  Fig.  18.  In  the  special  case  where  the 
misalignment  angle  0a  =  0°,  /3  is  an  integer,  and  there  is  no  intrinsic  clutter  motion, 
the  approximation  in  Eqn.  (100)  can  be  replaced  by  an  equality.  The  eigenspectrum 
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Pulse  0 
Pulse  1 
Pulse  2 

Element  Element  Element  Element 
0  12  3 

Figure  19  Effective  array  positions  for  successive  pulses.  Circles  indicate  indepen¬ 
dent  measurements. 

will  consist  of  exactly  Tc  nonzero  eigenvalues.  As  an  illustration  of  why  Brennan’s 
rule  holds,  consider  Fig.  19,  for  which  N  =  4,  M  =  3,  and  /3  =  1.  The  location 
of  the  second  element  during  the  first  pulse  is  the  same  as  the  location  of  the  first 
element  on  the  second  pulse.  Assuming  the  clutter  does  not  change  from  pulse  to 
pulse,  the  two  measurements  will  be  the  same.  Since  the  two  measurements  are 
not  linearly  independent,  the  rank  of  the  clutter  covariance  matrix  will  be  reduced. 
For  this  example  there  are  six  distinct  measurements,  which  is  equal  to  the  rank 
predicted  by  Brennan’s  rule.  Plots  of  the  eigenspectra  under  various  misalignment 
angles  are  shown  in  Fig.  20  The  amount  of  misalignment  has  little  effect  on  the  clut¬ 
ter  eigenspectrum.  Any  nonzero  misalignment  effectively  doubles  the  rank  of  the 
clutter  covariance  matrix.  This  is  as  expected,  since,  once  the  misalignment  angle 
deviates  from  zero,  there  are  two  distinct  clutter  ridges  present,  one  from  the  front 
lobe  and  another  from  the  back  lobe. 

3.6.3  Intrinsic  Clutter  Motion.  Intrinsic  clutter  motion  (ICM)  is  a  result 
of  variations  in  the  clutter  reflectivity  from  pulse  to  pulse.  The  variations  result  from 
movement  of  the  clutter  source,  such  as  trees  moving  in  the  wind  or  waves  on  the 
ocean.  Pulse-to-pulse  instabilities  of  the  radar  system  have  a  similar  effect  and  are 
included  as  part  of  the  ICM. 
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Figure  20  Clutter  eigenspectra  for  misalignment  angles  of  (f>a  =  0°,  45°, 

and  90°. 
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Prom  Eqn.  (84),  the  clutter  space  time  snapshot  was  given  by 

Nc  Nc 

Xc  =  £  =  52  QJfcbfc  ®  Sik- 

fc=l  k=l 

When  intrinsic  clutter  motion  is  present,  the  scalar  a*  must  be  replaced  with  an 
Mxl  vector  ak  =  [ccfc^o;  Oik,i]  •  •  • ;  Offc.M-i])  where  0;*;,^  is  the  random  amplitude  from 
the  kth.  clutter  patch  during  the  mth  PRI.  The  space-time  snapshot  is  now  given  by 

Nc 

Xc  =  S(«fcObfc)<8)afc.  (101) 

fe=i 

The  fluctuations  of  the  clutter  will  be  modeled  as  a  wide-sense-stationary  pro¬ 
cess  with  a  Gaussian  Doppler  spectrum.  The  temporal  autocorrelation  of  the  clutter 
fluctuations  will  then  be  given  by 

7c(m)  =  E  {aj+mO!*}  =  exp  (102) 

and  the  M  x  M  covariance  matrix  of  the  fluctuations  will  be  given  by 

Tk  =  E  {afcaf }  =  Toeplitz[7c(0);  7c(l); . . . ;  7c(M  -  1)],  (103) 


The  spectral  standard  deviation  Kc  is  generally  expressed  in  terms  of  a  velocity 
standard  deviation  cr^  as 
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(104) 


The  velocity  standard  deviation  is  a  measured  quantity  for  various  conditions  and 
types  of  clutter.  With  these  deflnitions,  the  clutter  covariance  matrix  can  be  written 
as 


(105) 
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3.1  Fully  Adaptive  Space-Time  Adaptive  Processing  (STAP) 

A  fully  adaptive  space-time  adaptive  processor  computes  and  applies  a  new 
adaptive  weight  for  every  antenna  element  and  every  pulse.  Thus,  the  weight  vector 
is  of  size  MN.  The  data  space-time  snapshot  at  a  particular  range  is  given  by 

X  =  oii^t  +  (106) 

where  Vt  is  the  target  space-time  steering  vector  and  Xu  the  interference-plus-noise 
component  of  the  data.  The  optimum  weight  vector  can  be  calculated  from  [5] 

w  =  cR;;Vt,  (107) 

where  c  is  a  constant  and  is  the  interference-plus-noise  covariance  matrix  as 
given  in  the  preceding  sections.  This  weight  vector  is  optimum  in  the  sense  that  it 
maximizes  the  signal-to-interference-plus-noise  ratio  (SINK)  and  the  probability  of 
detection  for  a  given  probability  of  false  alarm,  and,  with  the  proper  choice  of  c, 
minimizes  the  the  output  power  subject  to  a  gain  constraint  in  the  target  direction. 
[8] 

Once  the  weight  vector  is  computed,  the  adapted  pattern  can  be  determined 

from 

Pw(i?)^)  =  |w^v(t?,'n7)  .  (168) 

Figure  21  shows  a  sample  array  pattern  of  the  beamformer.  The  values  of  the 
parameters  used  are  given  in  Table  2. 


43 


Table  2  Parameters  used  to  generate  Fig.  21 


Number  of  elements 

N 

15 

Number  of  pulses/CPI 

M 

15 

Number  of  clutter  patches 

Nc 

360 

Transmitter  peak  power 

Pt 

200  kW 

Operating  frequency 

fo 

450  MHz 

Pulse  width 

200  fxs 

PRF 

fr 

300  Hz 

Antenna  transmit  gain 

Gt 

15  dB 

Element  gain 

9 

4dB 

Element  pattern 

<i>) 

Cosine 

Instantaneous  Bandwidth 

B 

4  MHz 

Noise  figure 

No 

3dB 

System  Losses 

Ls 

4dB 

Platform  height 

ha, 

9000  m 

Platform  velocity 

Va 

50  m/s 

Velocity  misalignment  angle 

00 

0° 

Target  azimuth 

<i>t 

10° 

Target  Doppler 

ft 

100  Hz 

Target  range 

Rt 

130  km 

Target  RCS 

2  m2 

Clutter  terrain  parameter 

7 

-3  dB 

Number  of  jammers 

J 

2 

Jammer  ERPD 

Sj 

1000  W/Hz 

Jammer  range 

Rj 

370  km 

Jammer  azimuth  angles 

-35°,  25° 

Jammer  elevation  angles 

li _ 

o 

O 

o 

O 

Doppler  Frequency  (Hz) 


-80  -60  -40  -20  0  20  40  60  80 

Azimuth  Angle  (deg) 

Figure  21  The  adapted  pattern  in  the  presence  of  a  target  at  10°,  100  Hz,  two 
jammers  at  25°  and  —35°,  and  clutter. 
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3.8  Signal-to-Interference-plus-Noise  Ratio 

A  common  measure  of  performance  is  the  output  signal-to-interference-plus- 
noise  ratio  (SINR).  The  output  signal  from  the  beamformer  is  given  by 

z  =  +  (109) 


where  the  first  term  is  due  to  the  target  signal  and  the  second  term  is  due  to  the 
undesired  signals.  Computing  the  power  in  each  signal  and  finding  the  ratio  give  the 
SINR, 


SINR  = 

w^RuW 


w“jn.uw 

If  the  optimum  weight  vector  is  used,  this  expression  reduces  to 


(110) 


SINRopt  =  ^tyfli~W 


(111) 


By  letting  Vt(‘Ci7)  =  b(n7)  0  a(i?f),  the  SINR  can  be  computed  as  a  function  of  target 
Doppler  from 

57iVR(n7)  =  ^tVt(ct7)-®^R“^vt(n7).  (112) 

Similarly,  the  SINR  can  be  calculated  as  a  function  of  target  azimuth  as 

SINR{^)  =  (113) 

where  ^t(i?)  depends  on  the  gain  of  the  antenna  when  steered  in  the  direction  of 
interest. 


3.9  Summary 

In  this  chapter,  we  have  presented  a  theoretical  development  of  the  signals 
involved  in  space-time  adaptive  processing.  For  a  linear  array  of  N  elements  which 
transmits  a  coherent  burst  of  M  pulses,  the  spatial  and  temporal  steering  vectors 
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are  defined  to  be 


a(i?)  = 
b(zi7)  = 

and  the  space-time  steering  vector  is  defined  as  the  Kronecker  product  of  these  two 
vectors, 

v(i?,ti7)  =  b(n7)  0  &{'&). 

The  target  space-time  snapshot  is  found  by  multiplying  the  target  space-time  steering 
vector  by  the  complex  amplitude  of  the  target  return, 


Xt  =  Oit^t  ®  at. 


The  space-time  snapshot  of  the  noise  internal  to  the  receiver  is  modeled  as  a  random 
vector,  and  the  covariance  matrix  of  the  noise  was  represented  as  a  MiV  0  MN 
identity  matrix.  The  jammer  signal  is  modeled  as  a  random  vector  ,Q'j,  of  length 
M.  The  space-time  snapshot  for  the  jammer  is  the  Kronecker  product  of  this  signal 
vector  with  the  jammer  steering  vector, 

Xj  =  OLj  0  Hj. 

The  jammer  covariance  matrix  is  given  by 

Rj  =  Im  ® 

where  Aj  is  a  matrix  of  the  spatial  steering  vectors  for  all  the  jammers,  and  Sj  is  a 
diagonal  matrix  of  the  jammer  powers.  The  clutter  space-time  snapshot  is  the  sum 
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of  the  space-time  snapshots  of  all  the  clutter  patches, 

Nc 

Xc  =  £  ®  afc. 

*=1 

The  clutter  covariance  matrix  is  given  by 

R=  V  S  V 

C  —  »  c'-'c  c> 

where  Vc  is  a  matrix  of  the  space-time  steering  vectors  to  each  clutter  patch  and  Sc 
is  a  diagonal  matrix  of  the  clutter  power  distribution. 

The  total  covariance  matrix  is  the  sum  of  the  covariance  matrices  of  the  noise, 
jammer,  and  clutter  signal.  The  optimum  weight  vector  is  found  by  inverting  the 
covariance  matrix  of  the  undesired  signals  and  multiplying  by  the  target  spatial 
steering  vector, 

w  =  R;^vt. 

The  adapted  angle-Doppler  pattern  of  the  antenna  is  found  from 

Py,{'d,w)  =  w^v(i?,tt7)  . 

The  SINK  is  used  as  a  measure  of  the  processor.  If  the  optimum  weight  vector  is 
used,  the  SINK  is  given  by 

SINR  =  (,V«RZW 

The  following  chapter  extends  these  results  to  include  the  effects  of  near-field  scat- 
terers. 
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IV.  Near-Field  Scatterers 


4.1  Introduction 

To  this  point  in  the  analysis  of  space-time  adaptive  processing,  we  have  as¬ 
sumed  the  array  is  in  free  space.  This  is  clearly  not  a  realistic  situation.  The  array 
must  be  mounted  on  a  platform  of  some  type.  Depending  on  the  position  of  the 
array,  parts  of  the  platform,  such  as  the  wing  or  the  tail,  may  be  in  the  field  of 
view  of  the  array.  The  signals  from  the  target,  jammers  and  clutter  will  refiect  off 
these  near-field  scatterers  (NFS)  into  the  array,  increasing  the  complexity  of  the 
total  signal  received  by  the  array. 

Yu,  Burnside,  and  Gilreath  [10]  developed  a  geometrical  theory  of  diffraction 
(GTD)  model  of  the  signals  which  will  enter  an  antenna  mounted  on  the  surface  of  an 
airframe.  Their  model  included  the  effects  of  refiections  off  the  airframe  waves  which 
strike  the  airframe  at  a  tangent  and  propagate  along  the  surface  of  the  airframe  to 
the  antenna.  Barile  et  al.  [3]  describe  the  effect  of  a  near-field  scatterer  on  a  two- 
element  displaced  phase  center  antenna  (DPCA).  Figure  22  illustrates  the  four  paths 
by  which  a  signal  transmitted  by  the  first  element  can  travel  to  the  target  and  back 
to  the  array.  Figure  22(a)  shows  the  direct  path  for  the  transmitted  and  received 
signals.  This  is  the  same  as  in  Chapter  III.  Figure  22(b)  shows  the  case  in  which  the 
transmitted  signal  follows  the  direct  path  to  the  target  and  refiects  off  the  NFS  on 
return  before  entering  the  array.  Figure  22(c)  shows  the  transmitted  signal  reflecting 
off  the  NFS  before  traveling  to  the  target  and  following  the  direct  path  to  return  to 
the  array.  Figure  22(d)  shows  the  double-bounce  case,  in  which  the  both  the  transmit 
and  received  signals  reflect  off  the  NFS.  Each  combination  of  antenna  element  and 
near-field  scatterer  will  produce  a  similar  set  of  signal  paths.  For  the  purpose  of 
this  analysis,  the  signals  which  scatter  off  the  NFS  more  than  once,  as  in  Figure 
22(d),  will  be  assumed  to  have  negligible  power  compared  to  the  direct  path  and 
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Figure  22  Four  possible  paths  for  a  signal  transmitted  from  the  first  element  to 
travel  to  the  target  and  return  to  the  array:  (a)  direct  path  to  and  from 
target,  (b)  Path  B,  direct  path  to  target,  scattering  path  on  return,  (c) 
Path  C,  scattering  path  to  target,  direct  path  on  return,  (d)  Path  D, 
scattering  path  to  and  from  target. 

single  bounce  paths.  In  the  case  of  multiple  NFS,  this  assumption  will  also  eliminate 
signals  which  bounce  off  two  or  more  different  scatterers. 

In  a  later  paper  [4],  Barile  et  al.  discuss  the  effect  of  a  near-field  scatterer  on  the 
signal  incident  on  an  array  in  the  presence  of  range  ambiguous  clutter.  Each  range 
ambiguity  will  contain  a  component  plane  wave  with  the  same  Doppler  frequency  as 
the  target.  When  this  component  scatters  off  the  NFS  into  the  array,  part  of  the 
scattered  wave  will  have  the  same  horizontal  wavenumber  and  Doppler  frequency 
as  the  target  return,  but  a  different  vertical  wavenumber.  A  linear  array  like  the 
one  discussed  in  Section  3.2,  which  is  only  adaptive  in  azimuth,  will  not  be  able  to 
distinguish  between  the  target  and  the  interference.  Fully  adaptive  STAP  with  a  two- 
dimensional  array  which  is  adaptive  in  both  azimuth  and  elevation  will  accurately 
null  out  the  signal  from  the  NFS,  but  at  significant  cost  in  processing  time.  In  [4], 
the  authors  propose  preadaptive  near-field  nulling  to  reduce  the  processing  required 
with  a  two-dimensional  array.  The  NFS  is  assumed  to  be  at  a  fixed,  measurable 
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Figure  23  Geometry  of  the  radar  system  including  the  near-field  scatterer. 

location  relative  to  the  array.  The  receiver  weights  required  to  place  a  null  in  the 
direction  of  the  NFS  are  precomputed  for  each  column  of  the  array.  The  adaptive 
processor  then  weights  the  rows  of  the  array  to  null  out  clutter  and  jamming  signals 
in  azimuth.  This  procedure  reduces  the  degrees  of  freedom  available  to  the  processor, 
but  also  reduces  the  complexity  of  the  calculations  required  to  shape  the  beam. 

This  chapter  will  parallel  the  derivation  presented  in  Chapter  3,  but  include 
the  effects  of  the  near-field  scatterer.  The  scatterers  will  be  assumed  to  be  point 
scatterers,  from  which  the  scattered  wave  pattern  will  be  isotropic  and  independent 
of  the  direction  of  incidence.  Because  the  scatterer  is  in  the  near-field  of  the  array, 
the  scattered  waves  must  be  treated  as  spherical,  meaning  each  receiver  in  the  array 
will  detect  a  signal  with  different  amplitude  and  direction  of  arrival. 
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4-2  Radar  System 

The  geometry  of  the  system  is  the  same  as  described  earlier,  with  the  addition 
of  the  near-field  scatterer  (NFS),  as  illustrated  in  Fig.  23.  To  begin,  a  single  near¬ 
field  scatterer  is  located  in  polar  coordinates  at  (Zog,  <j>os),  where  the  origin  is  the  0th 
element  of  the  array  and  the  angle  (j)  is  measured  from  the  normal  to  the  array.  The 
target  signal  originates  from  an  angle  <f)t,  and  a  jammer  signal  is  at  0^.  The  antenna 
elements  are  identical  and  are  equally  spaced  along  the  aj-axis.  The  NFS  is  assumed 
to  be  a  point  scatterer  close  enough  to  the  array  to  require  the  scattered  wave  to 
be  treated  as  spherical.  The  NFS  is  assumed  to  have  no  effect  on  the  temporal 
information  in  the  plane  wave,  so  the  temporal  steering  vector  b  is  the  same  for 
the  scattered  and  direct  path  waves.  Because  the  scatterers  are  close  to  the  array 
(within  a  few  meters),  the  delay  of  the  scattered  signal  is  be  small  compared  to  the 
pulse  width.  Therefore,  the  scattered  signal  enters  the  array  in  the  same  PRI  as  the 
direct  path  signal. 

The  distance  from  the  scatterer  to  the  nth  antenna  element  is  given  by 

Ins  =  [4  +  -  2losndsm{(f)ns)]^^^-  (114) 

For  the  case  of  multiple  scatterers,  the  subscript  s  denotes  the  scatterer  being  refer¬ 
enced,  so  that  Ins  is  the  distance  from  the  nth  antenna  element  to  the  sth  scatterer. 
The  projection  of  Iqs  onto  the  direct  path  to  the  target  is  given  by 


I'ts  =  hs  cos((l>t  -  (l>os)-  (115) 

For  the  direct  path  signal,  the  incident  wave  was  assumed  to  be  planar,  so  the 
signal  received  by  each  element  contained  the  same  amount  of  energy  and  the  phase 
progression  across  the  array  was  linear.  The  wave  from  a  near-field  scatterer  is  not 
planar.  Each  element  receives  a  different  amount  of  scattered  energy,  determined  by 
its  distance  from  the  scatterer.  The  ratio  of  the  power  of  the  scattered  signal  received 
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Table  3  Symbols  for  near-field  scatterer 


<t>t 

angle  to  target 

4*j 

angle  to  jammer 

i>d 

phase  shift  imposed  by  NFS 

^ns 

angle  from  nth  element  to  sth  scatterer 

s 

number  of  near-field  scatterers 

^ns 

distance  from  nth  element  to  sth  scatterer 

=  [/g,  -f-  {ndf  -  2lQsnds\rx{(j)ns)]^'^ 

I'ts 

projection  of  Iqs  onto  direct  path 

“  los 

(^S 

NFS  bistatic  radar  cross  section 

(ns 

scattering  loss 

by  the  nth  element  to  the  power  of  the  signal  incident  on  the  NFS  will  be  referred 
to  as  the  scattering  loss  and  is  given  by  (ns  =  ^sH^T^ln^i  where  Gs  is  the  bistatic 
radar  cross  section  (RCS)  of  the  scatterer.  For  a  point  scatterer,  Cg  is  independent 
of  incident  angle  and  the  scattered  wave  will  be  isotropic.  The  power  of  the  signal 
from  the  scatterer  incident  on  the  nth  antenna  element  is  found  by  multiplying  the 
power  incident  on  the  scatterer  by  the  scattering  loss. 

Because  the  distance  from  the  scatterer  to  each  element  in  the  antenna  array 
has  a  nonlinear  dependence  on  the  interelement  spacing,  the  phase  progression  of  the 
steering  vector  is  nonlinear.  The  relative  phase  shift  at  the  nth  element,  referenced 
to  the  0th  element,  will  be  given  by  (27r/Ao)(ios  ~  Ins)-  There  will  also  be  a  phase 
shift  of  =  180°  when  the  signal  reflects  off  the  NFS.  The  parameters  related  to 
the  near-fleld  scatterer  are  summarized  in  Table  3. 

The  space-time  snapshot  will  be  similar  in  form  to  Eqn.  (50), 


X  =  Xt  +  Xu, 


(116) 


where  the  tilde  (  ~  )  indicates  the  presence  of  signals  from  the  NFS.  In  addition  to 
the  direct  path  terms,  the  new  target  space-time  snapshot  will  include  terms  due  to 


53 


the  scattered  target  signals,  and  the  undesired  signal  component  will  include  terms 
due  to  the  scattered  clutter  and  jammer  signals. 

The  scattered  target  signal  does  not  have  the  same  spatial  steering  vector  as  the 
direct  path  signal,  but  the  temporal  steering  vector  is  unaffected  and  the  amplitude 
of  the  signal  is  still  be  deterministic.  Therefore,  the  scattered  target  signal  only  shifts 
the  mean  of  the  space-time  snapshot,  leaving  the  covariance  of  the  data  unaffected. 
The  total  covariance  matrix  is  now  given  by 

^  =  E[(XuXu] 

=  E[(Xc  +  Xj  +  Xn)iXc  +  Xj+Xnn  (117) 

The  direct  path  and  scattered  clutter  returns  are  correlated,  and  therefore  must  be 
combined  in  the  computation  of  the  clutter  covariance  matrix.  The  same  is  be  true 
of  the  jammer  signal.  The  clutter,  jammer,  and  noise  signals  are  assumed  to  be 
mutually  uncorrelated,  so  the  total  covariance  matrix  can  be  written  as 

^  =  E[XcXc]  +  ^IXjXf]  +  E[XnXn] 

—  R,c  R.J  -l-  Rji,  (11®) 

where  the  tilde  (  ~  )  indicates  the  covariance  matrix  includes  the  effects  of  signals 
from  the  near-field  scatterer. 

4.3  Scattered  Target  Signal 

The  signal  received  from  the  target  by  the  nth  antenna  element  is  the  sum  of 
the  direct  and  scattered  paths.  For  the  mth  pulse,  the  direct  path  signal  at  the  nth 
element  is  given  by  Eqn.  (64)  and  is  repeated  here: 
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The  signal  which  travels  the  direct  path  to  the  target  and  bounces  off  the  NFS  on 
the  return  (Path  B)  is  indicated  by  the  minus  (-)  superscript  and  is  given  by 

^mn  =  (119) 

where  is  as  defined  in  Eqn.  (115)  and 

1 

(120) 

This  expression  is  obtained  by  computing  the  single  pulse  signal-to-noise  ratio  (SNR) 
for  the  scattered  signals, 

(ATrYNoLsRt  ^  ^ 

and  dividing  the  result  by  the  SNR  of  the  direct  path  signal  given  by  Eqn.  (65).  All 
the  terms  will  cancel  except  the  element  gains  in  the  target  and  scatterer  directions, 
g(<(>t)  and  g{(f)ns)i  and  the  scattering  loss,  Cns-  Equation  (120)  represents  the  energy 
ratio  of  the  scattered  path  signal  to  the  direct  path  signal.  The  subscripts  t,  n  and 
s  represent  respectively  the  target,  the  nth  antenna  element,  and  the  sth  near-field 
scatterer.  By  defining 

<«<;.=  (122) 

Eqn.  (119)  can  be  written  as 

(123) 

This  expression  is  similar  to  Eqn.  (64),  with  replacing  The  space- time 

snapshot  of  the  signals  which  have  followed  the  path  of  Fig.  22(b)  can  be  written  as 

Xt"  =  afbt  ®  dt",  (124) 


^tns 


g{<t>\ 


I  9{<l>t) 
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where  •  •  • ;  represents  the  spatial  steering  vector  for 

the  scattered  target  signal.  Because  each  element  is  a  different  distance  from  the 
scatterer,  each  element  of  this  vector  has  a  different  amplitude,  and  the  phase  pro¬ 
gression  across  the  elements  is  nonlinear. 

The  signal  entering  the  nth  receiver  which  bounced  off  the  NFS  before  going 
to  the  target  [Fig.  22(c)]  is  indicated  by  the  plus  (-f-)  superscript.  This  signal  is 
more  complicated.  Because  each  element  is  a  different  distance  from  the  NFS,  there 
will  appear  to  be  N  pulses  being  transmitted  from  the  NFS,  each  with  a  different 
amplitude  and  phase.  The  signal  received  by  the  nth  element  is  the  sum  of  these  N 
pulses: 


Xtn  =  0‘t 


HU 


-l-n+  ^ 

N-1 


^m2irwi  gjfn27riJ<  gii&d 


p=0 

JV-1 

=  X)  <^tp$ 

p=0 


(125) 


where 


^tps 

9{^ps)  >  ^ 

(126) 

^tps 

(127) 

N-l 

^tns 

_  gjn27n2tgj^d  ^ 

(128) 

p=0 


The  space-time  snapshot  due  to  the  signals  in  Eqn.  (125)  can  be  written  as 


Xt  =  tttbt  0  dt+. 


(129) 
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The  total  target  space-time  snapshot  is  found  by  adding  Eqns.  (124),  (129), 


and  (71) 


Xt  =  Xt  +  Xt  +Xt 

=  oit{ht  ®  at)  +  Q!t(bt  <2)  dj“)  -I-  Q!t(bt  <S)  d^) 

=  Q!fbt  0  (at  -t-  d^j  -I-  d^).  (130) 


For  multiple  scatterers,  Eqns.  (122)  and  (128)  must  be  summed  over  all  scatterers, 
and  Eqn.  (130)  becomes 


Xt 


atbt  0 


+  52(dts  +  d^) 


s=0 


atbt  0  if 


(131) 


The  term  at  in  Eqn.  (131)  can  be  considered  the  spatial  steering  vector  of  the  target 
signal  entering  the  array  through  both  the  direct  and  scattered  paths.  The  Kronecker 
product  of  the  temporal  steering  vector,  bt,  with  this  spatial  steering  vector  gives 
the  space-time  steering  vector  for  the  target  due  to  both  the  direct  and  scattered 
paths, 

Vt  =  bt  0  at.  (132) 

This  expression  will  be  used  later  in  the  calculation  of  the  ideal  weight  vector. 


4-4  Scattered  Jammer  Signal 

We  will  begin  by  considering  a  single  jammer  and  single  NFS.  The  parameter 
Ijg  is  similar  to  given  by  Eqn.  (115),  except  the  angle  to  the  jammer  is  used  in 
place  of  the  target  angle, 

Ijg  =  Zqs  cos(^j  00s)"  (^33) 

The  total  jammer  signal  contains  two  components:  the  direct  path  and  the  Jammer- 
NFS-array  path.  During  the  mth  PRI,  the  signal  received  by  the  nth  receiver  which 
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has  scattered  off  the  NFS  will  be  given  by 


ran  '^rnf*jns^  ^  > 


(134) 


where 

(135) 

Prom  these  expressions,  the  components  of  the  spatial  steering  vector  of  the  scattered 
jamming  signal  can  be  written  as 


^jns  — 


9i^ns) , 


ilr,.  =  a- (136) 

By  defining 

kj  =  aj  +  (137) 

the  jammer  space-time  snapshot  is  now  given  by 

Xj  =  aj  0  aj  =  OLj  0  -t-  otj  0  =  Xj  +  Xjs'  (138) 


The  jammer  covariance  matrix  is  computed  as  in  Eqn.  (78),  with  kj  replacing  a^. 


E  {(ttj  0  kj){aj  0  kj)^} 
E{ajaf)  0  (a^af ) 

®  (aj-af ) 

1m  <0 


(139) 


For  J  jammers,  can  be  extended  in  the  same  manner  as  in  Section  3.5, 


A  W  - A-®^ 


(140) 
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where 


=  [ai,  §2,  •••,  a^].  (141) 

For  S  scatterers,  djg  in  Eqn.  (137)  must  be  summed  over  all  the  scatterers, 

aj  =  aj  +  E)  (142) 

3=0 

This  expression  can  then  be  used  in  Eqn.  (141). 

4-5  Scattered  Clutter  Signal 

The  return  from  each  clutter  patch  consists  of  a  set  of  signals  similar  to  the  set 
of  target  signals.  There  is  the  direct  path  signal,  the  two  single-bounce  path  signals, 
and  a  double-bounce  path.  As  with  the  target  signal,  the  double-bounce  path  will  be 
ignored  in  computing  the  total  clutter  return.  It  is  assumed  there  is  no  misalignment 
angle  {(pa  =  0)  and  no  intrinsic  clutter  motion. 

The  direct  path  signal  is,  from  Eqn.  (84), 

^e  Nc 

Xc  =  E  a'fcVjfc  =  2  «*(b*  ®  a*). 
k=l  fc=l 

The  signal  which  bounces  off  the  NFS  while  returning  from  the  feth  clutter  patch 
[Fig.  22(b)]  is  given  by 

=  (143) 


59 


The  space-time  snapshot  for  this  set  of  signals  is  found  by  summing  over  all  the 
near-field  scatterers  and  all  the  clutter  patches, 

Nc  /  S-1  \ 

Xj  =  (144) 

k=l  \  s=0  / 

The  signal  which  bounces  off  the  NFS  before  traveling  to  the  fcth  clutter  patch 
[Fig.  22(c)]  is  given  by 


X 


+  _ 


kmn 


=  Oik 


gjm27rn7t  g  jn27rtf  t  gj  V'd 


jm2nwt ^jn2ir'&k ^ipd  ^  q+ 


N-1 

p=0 

N-1 

p=0 


(145) 


The  space-time  snapshot  due  to  these  signals  is  given  by 

Xc=^0!k  (bk  0  £  dl")  .  (146) 

fc=l  \  s=0  / 

The  total  clutter  space-time  snapshot  is  found  by  adding  Eqns.  (144)  and  (146)  to 
(84), 


Xc  = 


.+ 


Xc  +  Xc  +Xc 

Nc  r 


S-l 


=  E 

A:=l  L 
Nc  / 

=  E  (  °‘k^k  ® 

k=l  \ 


Okbk  0  afc  -h  52  (ajfcbfc  ®  -h  afcbfc  0  djj,) 

5=0 

a*  +  E  (‘^fcs  +  ‘^fcs) 


s=0 


(147) 


By  defining 


y 


Figure  24  Location  of  the  near-field  scatterers  relative  to  the  array. 

and 

Vc=[vi;  V2;  •••;  v^J,  (149) 

the  total  clutter  covariance  matrix  can  be  written  in  the  same  form  as  Eqn.  (89), 

(150) 

The  clutter  covariance  matrix  will  not  be  Toeplitz-block-Toeplitz  as  in  Section  3.6. 
The  blocks  within  the  matrix  lose  their  Toeplitz  structure  when  the  near-field  scat¬ 
terers  are  introduced.  Since  the  temporal  information  of  the  signal  is  not  affected, 
the  matrix  will  retain  its  block-Toeplitz  structure.  The  clutter  covariance  matrix 
shown  in  Fig.  25  was  generated  for  the  same  conditions  as  used  in  Fig.  14,  but  with 
the  addition  of  seven  near-field  scatterers  with  a  RCS  of  0.5  m^,  spaced  equally,  one 
meter  in  front  of  the  array,  as  illustrated  in  Fig.  24.  Figure  26  is  a  magnified  view 
of  the  upper  left  corner  of  the  covariance  matrix 

4.6  Fully  Adaptive  Processing  in  the  Presence  of  Near- field  Scatterers 

As  in  Chapter  3,  the  total  covariance  matrix  is  the  sum  of  the  jammer,  clutter, 
and  noise  covariance  matrices.  In  the  presence  of  near-field  scatterers,  Eqns.  (139), 


Re  =  VeSeVf . 
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Figure  25  Clutter  covariance  matrix  in  the  presence  of  seven  near-field  scatterers. 
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Figure  26  Upper  left  corner  of  the  clutter  covariance  matrix  in  the  presence  of 
seven  near-field  scatterers. 
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(150),  must  be  used  in  place  of  Eqns.  (78)  and  (89), 


=  Rj  +  Rc  +  Rji- 


(151) 


The  form  of  the  optimum  weight  vector  is  the  same  as  Eqn.  (107),  but  using  the  co- 
variance  matrix  of  Eqn.  (151)  and  the  target  space-time  steering  vector  of  Eqn.  (132). 
The  optimum  weight  vector  is  given  by 


w  = 

w  = 


(152) 


where  and  are  as  defined  in  Eqns.  4.3.  This  weight  vector  is  then  used  to 
determine  the  adapted  angle- Doppler  pattern  of  the  antenna  as  in  Eqn.  (108), 


Pw(i?, 


w^v(i?,  w)  ^ . 


(153) 


The  resulting  adapted  pattern  is  shown  in  Fig.  27.  The  same  parameters  were  used 
as  listed  in  Table  2,  with  the  scatterer  configuration  shown  in  Fig.  24.  Note  how  the 
sidelobe  structure  has  been  distorted  in  the  presence  of  near-field  scatterers.  This 
distortion  causes  the  clutter  and  jammer  nulls  to  become  blurred,  allowing  more 
energy  from  these  sources  to  enter  the  receiver. 


4-7  Summary 

In  this  chapter,  the  theoretical  development  of  space-time  adaptive  processing 
presented  in  Chapter  3  has  been  extended  to  include  the  effects  of  near-field  scat¬ 
terers.  The  near-field  scatterers  have  no  effect  on  the  temporal  steering  vectors  of 
the  received  signals;  only  the  spatial  steering  vectors  are  affected.  While  the  spatial 
steering  vectors  for  the  direct  path  of  each  signal  have  unit  amplitude  and  a  lin¬ 
ear  phase  shift,  the  spatial  steering  vectors  of  the  scattered  signals  have  neither  of 
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Doppler  Frequency  (Hz) 


Azimuth  Angle  (deg) 

Figure  27  The  adapted  pattern  for  the  same  conditions  as  Fig.  21,  with  the  addi¬ 
tion  of  seven  near-field  scatterers,  as  illustrated  in  Fig.  24. 
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Relative  Power  (dB) 


these  qualities.  Because  each  scatterer  is  a  different  distance  from  each  element,  the 
signals  received  by  each  element  from  a  particular  scatterer  have  a  nonlinear  phase 
progression.  Also,  since  the  scatterers  are  in  the  near-field  of  the  array,  the  scattered 
wave  must  be  treated  as  spherical,  meaning  the  amplitude  of  the  scattered  signals 
received  by  each  element  is  different. 

To  simplify  the  computations,  signals  which  scatter  off  more  than  one  near¬ 
field  scatterer  are  ignored.  For  the  path  shown  in  Fig.  22(b),  in  which  the  signal 
bounces  off  the  scatterer  when  returning  to  the  array,  the  spatial  steering  vector  is 
given  by 

~  ■"> 


where 


“tn«  ~  “tns®  °  ^  > 


atns  is  the  ratio  of  the  amplitude  of  the  direct  path  target  signal  to  the  scattered 
path  target  signal  from  the  sth  scatterer  to  the  nth  element,  is  the  distance  from 
the  nth  element  to  the  sth  scatterer,  and  is  the  length  of  the  projection  of 
onto  the  target  vector.  The  spatial  steering  vector  of  the  signal  which  bounces  off 
the  scatterer  before  going  to  the  target  is  given  by 


[^tOs  > 


f]+ 


where 

p=0 

The  total  spatial  steering  vector  for  the  target  signal  is  the  sum  of  the  spatial  steering 
vectors  for  the  direct  path  and  scattered  paths  for  all  the  scatterers 


( 


s-i 

at  =  at  -t-  53 


s=0 
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this  expression  is  then  used  to  compute  the  space-time  steering  vector  and  space-time 
snapshot  of  the  target  signal  as  in  Chapter  3, 


Vi  =  bt  (g)  at. 


The  scattered  jammer  signal  has  the  same  form  as  the  target  signal  which 
bounced  oflP  the  scatterer  while  returning  from  the  target.  The  components  of  the 
spatial  steering  vector  of  the  scattered  jammer  signal  are  given  by 


^jna 


where  and  have  the  same  form  as  for  the  target  signal,  but  using  the  direction 
to  the  jammer.  The  total  jammer  steering  vector  is  found  by  summing  dj^  over  all 
the  scatterers  and  adding  the  direct  path  steering  vector, 

s-i 

~ 

s=o 

The  jammer  covariance  matrix  is  computed  in  the  same  manner  as  in  Chapter  3, 
but  using  Sij  in  place  of  aj, 


where  Aj  =  [ai,  §2,  •  •  • ,  aj]. 

The  steering  vectors  for  the  scattered  clutter  returns  from  the  fcth  clutter  patch 
have  the  same  form  as  the  target  return. 


^kns 


a 


kna^ 


N-l 


p=0 
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The  total  spatial  steering  vector  for  the  kth  clutter  patch  is  the  sum  of  the  direct 
path  steering  vector  with  the  sum  of  the  scattered  path  steering  vectors  for  all  the 
scatterers, 

s-i 


ajb  =  ajb  +  5;  (dfc,  +  dl)  . 
8—0 


The  total  spatial  steering  vector  for  the  A:th  clutter  patch  becomes 


and  the  clutter  covariance  can  be  computed  as 


where  Vc  =  [vi,  V2,  •  • ' ,  ^Ncl- 

The  total  covariance  matrix  of  the  undesired  signals  is  the  sum  of  these  two 
covariance  matrices  with  the  noise  covariance  matrix, 

=  Rj  +  Rc  +  Rn* 

The  optimum  weight  vector,  adapted  antenna  pattern,  and  SINK  can  be  computed 
using  these  quantities, 


w 

SINK 


= ,  R;^v, 
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V.  Results 


5.1  Test  Plan 

To  determine  the  effects  of  the  near-field  scatterers  on  space-time  adaptive 
processing,  several  test  cases  were  examined.  The  number  of  scatterers  varied  from 
two  to  fifteen  and  the  RCS  of  the  scatterers  varied  from  0.01  to  0.5  m^.  Six  different 
orientations  were  used,  which  are  shown  in  Fig.  28.  For  orientations  A  through  E, 
the  scatterers  are  linearly  spaced  between  the  endpoints.  Orientation  F  is  similar  to 
orientation  A,  but  the  positions  of  the  scatterers  are  allowed  to  vary  by  ±  0.1  m  in 
both  the  X  and  y  directions.  Table  4  describes  the  orientations  in  more  detail.  In  all 
cases,  the  positions  of  the  scatterers  are  with  respect  to  the  left  end  of  the  array. 

5.2  Results 

5.2.1  No  External  Interference.  For  the  first  three  sets  of  figures,  the 
jammer  and  clutter  signals  were  ignored.  In  the  absence  of  external  interference,  the 
signal-to-noise  ratio  (SNR)  does  not  depend  on  the  Doppler  frequency. 

Figures  29  through  32  show  how  changing  the  RCS  of  the  scatterers  affects 
the  SNR.  The  scatterers  were  arranged  in  orientation  A.  Figure  29  shows  the  SNR 
as  a  function  of  azimuth  angle  and  scatterer  RCS  for  five  near-field  scatterers.  The 


Table  4 

Definition  of  ori 

entations 

Location  of 

Location  of 

Orientation 

first  scatterer 

second  scatterer 

(xi,yi)(m) 

(a;2,y2)(ni) 

A 

(0.1) 

(IM) 

B 

(1/2, 1) 

(1/2,1 +  1) 

C 

(0,1) 

(0,1) 

D 

(0,1) 

(1,1  +  1) 

E 

(0,1) 

(1, 1/2  +  1) 

F 

(0  ±  0.1, 1  ±  0.1) 

(1±0.1,1  +  1±0.1) 

is  the  length  of  the  array 


69 


: 

t 

I 

I 


(C)  (D) 


(E)  (F) 


Figure  28  Near-field  scatterer  orientations. 

SNR  decreases  as  the  scatterer  RCS  increases,  but  tends  to  fiatten  out  as  more 
scatterers  are  added.  For  larger  numbers  of  scatterers,  the  SNR  increases  as  the 
size  of  the  scatterers  increases.  The  increase  in  SNR  is  due  to  the  larger  scatterers 
refiecting  more  target  energy  into  the  array.  Figures  30  through  32  show  the  SNR  as 
a  function  of  azimuth  angle  and  scatterer  RCS  for  seven,  eight,  and  fifteen  scatterers 
respectively. 

The  notch  at  0°  in  Fig.  31  is  due  to  the  fact  that  each  NFS  is  directly  in  front 
of  an  antenna  element.  Since  the  NFS  imposes  a  180°  phase  shift  on  the  signal,  the 
signal  from  a  scatterer  to  the  element  directly  behind  it  will  be  180°  out  of  phase 
with  the  direct  path  signal.  The  sum  of  this  scattered  signal  with  the  direct  path 
signal  reduces  the  apparent  amplitude  of  the  direct  path  signal  entering  the  element. 
Because  only  eight  of  the  fifteen  elements  are  affected,  a  target  at  0°  appears  to 
produce  a  wavefront  with  an  amplitude  which  varies  with  position.  This  effect  is 
shown  in  Fig.  33.  This  problem  does  not  occur  with  15  scatterers  because  the  signal 
at  each  element  is  reduced  by  the  same  amount. 
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Figure  29  SNR  as  a  function  of  azimuth  angle  and  scatterer  RCS  for  5  NFS,  ori¬ 
entation  A. 
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Figure  30  SNR  as  a  function  of  azimuth  angle  and  scatterer  RCS  for  7  NFS,  ori¬ 
entation  A. 
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Figure  31  SNR  as  a  function  of  azimuth  angle  and  scatterer  RCS  for  8  NFS,  ori¬ 
entation  A. 
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Figure  33  Effect  of  eight  NFS  on  the  direct-path  wavefront  perceived  by  a  fifteen 
element  array. 


Figures  34  through  39  show  how  changing  the  number  and  orientation  of  the 
scatterers  affects  the  SNR.  For  each  orientation  in  Table  4,  the  target  angle  was 
varied  from  —50°  to  4-50°  and  the  number  of  scatterers  was  varied  from  zero  to 
fifteen.  The  scatterer  RCS  was  held  constant  at  0.5  m^.  For  orientation  A,  shown  in 
Fig.  34,  the  SNR  decreases  until  the  number  of  scatterers,  S,  reaches  five,  at  which 
point  the  SNR  begins  to  increase.  At  S  —  10,  the  SNR  is  at  the  same  level  as  for 
no  scatterers.  The  SNR  for  orientation  B  is  shown  in  Fig.  35.  The  minimum  SNR 
occurs  at  S  =  7,  and  the  SNR  returns  to  the  same  level  as  with  no  scatterers  at 
5  =  13. 

Figure  36  shows  the  SNR  for  orientation  C.  For  target  angles  greater  than  20°, 
the  scatterers  have  almost  no  effect  on  the  SNR  because  the  scattered  waves  do  not 
enter  the  array  from  the  same  direction  as  the  direct  path  wave.  For  targets  to 
the  left  of  the  array,  the  scatterers  affect  the  SNR  in  a  manner  similar  to  that  of 
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orientation  B.  The  SNR  decreases  as  the  number  of  scatterers  increases  from  zero  to 
seven,  then  increases  as  the  the  number  of  scatterers  continues  to  increase. 

Figures  37  and  38  show  the  SNR  plots  for  orientations  D  and  E.  As  the  angle 
between  the  scatterers  and  the  array  decreases,  the  SNR  plot  more  closely  resembles 
that  of  Fig.  34  for  orientation  A.  Figure  39  shows  the  SNR  for  orientation  F.  The 
plot  is  very  similar  to  Fig.  34,  indicating  the  precise  spacing  of  the  scatterers  has 
less  effect  on  the  SNR  than  their  number  and  general  configuration. 

Figures  40  through  42  show  the  adapted  pattern  as  a  function  of  azimuth 
angle  at  the  target  Doppler  for  different  numbers  of  scatterers  and  a  scatterer  RCS 
of  0.5  m^.  Figure  40  shows  the  adapted  azimuth  antenna  patterns  for  a  target  at 
0°,  100  Hz,  and  various  numbers  of  scatterers.  Fig.  41  is  for  a  target  at  20°,  and 
Fig.  42  is  for  a  target  at  40°.  The  shape  of  the  mainbeam  and  first  two  side  lobes 
is  not  significantly  affected  by  the  scatterers.  The  remaining  side  lobes,  however, 
show  some  significant  distortion,  especially  for  the  cases  of  five,  seven,  eight,  and  ten 
scatterers.  The  scatterers  are  spaced  more  than  A/2  apart,  and  the  criterion  to  avoid 
spatial  frequency  aliasing  given  in  Eqn.  (20)  is  not  met.  The  large  grating  lobes,  such 
as  in  Fig.  42  at  —20°  for  eight  scatterers,  are  the  result  of  this  spatial  undersampling. 
For  S  =  15,  the  side  lobes  are  lower  than  for  the  case  of  no  scatterers. 
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SNR  (dB) 


Figure  34  Signal-to-noise  ratio  as  a  function  of  azimuth  angle  and  number  of  scat- 
terers  using  orientation  A  and  tr,  =  0.5  m^. 


77 


SNR  (dB) 


50 


Number  of  scatterers 


0  -50 


Azimuth  Angle  (deg) 


Figure  35  Signal-to-noise  ratio  as  a  function  of  azimuth  angle  and  number  of  scat¬ 
terers  using  orientation  B  and  o-j  =  0.5  m^. 
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Figure  36  Signal-to-noise  ratio  as  a  function  of  azimuth  angle  and  number  of  scat¬ 
terers  using  orientation  C  and  £7,  =  0.5  m^. 
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SNR  (dB) 


Figure  37  Signal-to-noise  ratio  as  a  function  of  azimuth  angle  and  number  of  scat- 
terers  using  orientation  D  and  as  —  0.5  m^. 
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Relative  power  (dB) 


Relative  power  (dB) 


Relative  power  (dB) 


Table  5  Rank  of  the  clutter  covariance  matrix  for 
each  orientation  and  number  of  scatterers  S. 


5 

Scatterer  Orientation 

A 

B 

C 

D 

E 

F 

2 

82 

83 

83 

84 

84 

3 

83 

90 

90 

92 

90 

4 

86 

94 

93 

100 

95 

91 

5 

87 

96 

96 

105 

99 

93 

6 

87 

96 

97 

108 

101 

95 

7 

86 

95 

97 

109 

103 

95 

8 

84 

94 

97 

108 

102 

96 

9 

87 

95 

98 

111 

103 

97 

86 

95 

97 

113 

104 

98 

11 

85 

95 

96 

112 

104 

98 

12 

85 

94 

95 

110 

104 

97 

13 

85 

94 

95 

110 

104 

97 

14 

82 

93 

95 

110 

103 

98 

15 

79 

90 

94 

108 

103 

96 

5.2.2  External  Interference  Present.  The  following  plots  show  the  effects  of 
the  near-field  scatterers  on  the  adapted  antenna  pattern,  signal-to-interference-plus- 
noise  ratio  (SINK),  and  the  rank  and  eigenvalues  of  the  clutter  covariance  matrix. 
The  radar  system,  jammer,  and  clutter  parameters  used  are  the  same  as  those  listed 
in  Table  2. 

Table  5  shows  the  rank  of  the  clutter  covariance  matrix  for  a  target  at  0° 
and  each  combination  of  orientation  and  number  of  scatterers.  The  rank  indicates 
the  number  of  degrees  of  freedom  required  to  null  out  the  clutter  return.  The  two 
diagonal  orientations,  D  and  E,  have  the  greatest  effect  on  the  rank  of  the  clutter 
covariance  matrix.  Figure  43  shows  the  eigenvalues  for  each  orientation  with  eight 
scatterers.  In  most  cases,  the  highest  rank  occurred  with  between  seven  and  ten 
scatterers.  For  comparison,  the  rank  of  the  clutter  covariance  matrix  with  no  near¬ 
field  scatterers  is  29.  Increasing  the  scatterer  RCS  or  the  target  angle  had  little  effect 
on  the  rank. 
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Eigenvalue  number 


Figure  43  Eigenvalues  of  the  clutter  covariance  matrix  for  different  orientations  of 
8  scatterers,  target  at  0°,  ctj  =  0.5  m^.  The  first  solid  line  from  the 
left  represents  the  eigenvalues  of  the  clutter  covariance  matrix  with  no 
scatterers,  the  second  solid  line  is  for  orientation  A,  and  the  last  solid 
line  is  for  orientation  D. 
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Figures  44  through.  49  show  the  azimuth  antenna  pattern  at  the  target  Doppler 
for  different  orientations  and  numbers  of  scatterers.  The  target  is  at  20°,  100  Hz,  and 
the  scatterer  RCS  is  cxs  =  0.5  m^.  Appendix  B  contains  the  complete  angle-Doppler 
patterns  from  which  these  figures  were  taken.  The  presence  of  clutter  causes  the 
SINK  to  be  a  function  of  Doppler  frequency,  as  well  as  azimuth  angle.  Figures  50 
through  55  show  the  corresponding  SINK  vs.  Doppler  plots.  The  large  null  at  50 
Hz  corresponds  to  the  Doppler  frequency  of  clutter  at  20°.  The  smaller  nulls,  such 
as  the  one  at  -80  Hz  for  S=5  in  Fig.  50,  correspond  to  the  Doppler  frequencies  of 
clutter  entering  the  receiver  through  the  grating  lobes. 
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F.  The  number  of  scatterers,  S,  is  indicated  at  the 
The  dashed  line  is  the  SINK  with  no  scatterers. 


VI.  Conclusions 


6.1  Conclusions 

This  thesis  has  presented  a  derivation  of  space-time  adaptive  processing  in  the 
presence  of  near-field  scatterers.  In  Chapter  4,  we  derived  an  expression  for  the 
spatial  steering  vector  of  the  signal  entering  the  array  due  to  a  particular  scatterer. 
Summing  this  result  over  all  the  scattering  paths  and  adding  the  direct  path  steering 
vector  gives  the  steering  vector  for  the  combined  direct  and  scattered  signals.  This 
total  steering  vector  replaces  the  direct  path  steering  vector  in  the  computation  of 
the  space-time  snapshot  of  the  signal,  the  covariance  matrix  of  the  undesired  signals, 
the  optimum  weight  vector,  and  the  adapted  pattern  and  SINK. 

The  effect  of  near-field  scatterers  on  the  performance  of  space-time  adaptive 
processing  depends  on  the  number  and  configuration  size  of  the  scatterers.  The 
scatterers  increase  the  energy  entering  the  receivers  from  each  signal  source.  In 
some  cases,  the  fully  adaptive  processor  is  able  to  suppress  the  scattered  jammer 
and  clutter  energy  while  using  the  scattered  target  energy  to  increase  the  SINK. 
When  the  number  of  scatterers  is  small  and  the  spacing  between  them  is  more  than 
a  half  wavelength,  the  scatterers  create  grating  lobes  in  the  adapted  pattern.  The 
grating  lobes  allow  undesired  signals  to  enter  the  receiver,  reducing  the  SINK  by  as 
much  as  6  dB  in  some  situations.  As  the  number  of  scatterers  increases,  the  grating 
lobes  vanish. 

Changing  the  orientation  of  the  scatterers  had  the  most  dramatic  effect  on 
the  shape  of  the  adapted  antenna  pattern.  Arranging  the  scatterers  normal  to  the 
array,  as  in  orientations  B  and  C,  caused  the  greatest  disruption  of  the  sidelobe 
pattern,  in  some  cases  elevating  the  sidelobes  by  10  dB  or  more,  as  seen  in  Figs.  45 
and  46.  Although  the  elevated  sidelobes  allow  additional  undesired  signals  to  enter 
the  receiver,  the  adaptive  processor  is  able  to  maintain  good  SINK  performance. 
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experiencing  only  a  2-3  dB  loss  for  orientation  B  (Fig.  51)  and  a  achieving  a  slight 
gain  for  orientation  C  (Fig.  52). 

The  rank  of  the  clutter  covariance  matrix  when  the  near-field  scatterers  are 
present  is  significantly  higher  than  without  the  scatterers.  In  Section  3.6.2,  the 
return  from  a  clutter  patch  was  assumed  to  remain  constant  from  pulse  to  pulse. 
With  the  near-field  scatterers,  this  assumption  does  not  hold.  The  scattered  signals 
entering  one  element  will  be  diff'erent  from  the  scattered  signals  entering  the  next 
element  when  the  array  is  shifted,  as  in  Fig.  19.  The  maximum  values  for  the  rank 
for  each  orientation  occurs  with  between  eight  and  twelve  scatterers.  Increasing  the 
number  of  scatterers  beyond  5  has  relatively  little  effect  on  the  rank  when  compared 
to  the  increase  in  rank  when  the  first  two  scatterers  are  introduced. 

6.2  Suggestions  for  Future  Research 

In  every  case  considered  in  this  thesis,  the  target,  jammer,  and  scatterers  were 
all  in  the  same  plane.  The  effect  of  placing  the  scatterers  in  different  arrangements 
in  three  dimensions  should  be  studied.  The  extension  of  the  derivation  presented 
here  is  straightforward,  requiring  only  the  antenna  pattern  and  spatial  frequency  to 
be  calculated  as  a  function  of  elevation. 

The  effects  of  signals  bouncing  off  multiple  scatterers  also  needs  to  be  in¬ 
vestigated.  Multiple  scattering  will  probably  have  the  greatest  impact  at  angles 
approaching  ±  90°,  since  the  direct  path  gain  is  nearly  zero  at  these  angles. 
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Appendix  A.  Notational  Conventions  and  Matrix  Operations 

This  Appendix  defines  the  notation  conventions  which  will  be  used  in  this  paper. 

Vectors  (defined  to  be  column  vectors):  lower-case  bold  letters:  x,  a. 

When  the  elements  of  a  vector  are  given  explicitly,  a  semicolon  (;)  between 
elements  indicates  a  column  vector,  and  a  comma  (,)  or  no  punctuation  between 
elements  indicates  a  row  vector. 

Row  Vector  b  =  [ci,  02,  ...,  On] 

=  [fli  0,2  ...  Un] 

Column  Vector  a  =  [ai;  02;  . . . ;  ««] 

=  a 

Matricies:  upper-case  bold  letters:  R,  S. 

Element  of  matrix  A  =  ak,i 
Scalars:  unbold  characters:  a,  N. 

Complex  conjugate:  (•)* 

Transform  of  a  vector  or  matrix:  (•)^ 

Conjugate  Transform:  (•)^ 

M  X  M  Identity  matrix: 

M  xl  vector  of  ones:  lAf=[l;  1;  •  •  • ;  1] 

Toeplitz  matrix:  a  matrix  in  which  all  diagonal  elements  are  equal.  The  vectors 
X  and  y  represent  the  first  row  and  first  column  respectively.  (In  this  example. 
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M  =  N,  but  this  is  not  necessary.) 


/ 

\ 

Xi 

X2 

X3 

Xn 

2/2 

Xi 

X2  ••• 

Xn-1 

Toeplitz(x,  y)  = 

2/3 

y2 

Xi  ••• 

Xn-2 

^  Vm 

2/m-i 

yM-2  •  •  • 

Xl  j 

Toeplitz(x)  the  same  vector  is  used  for  the  first  row  and  first  column,  and 
therefore  the  resulting  matrix  is  symmetric. 

Inner  product  of  x  and  y,  each  of  length  M:  z  =  x^y  =  x  •  y  =  aibi  +  02&2  + 

- 1- 

Outer  product  of  x  and  y  of  lengths  M  and  N  respectively: 


Z  =  xy^  = 


fll6i  O162  *  *  ■ 

0261  01262  ’  ’  ■  o,2bi^ 


(lM^2  '  •  * 


Kronecker  Product:  for  A  (M  x  N)  and  B  (if  x  L).  The  resulting  matrix  Z  is 
MK  X  NL. 

r 

oi.iB  aifi'B  •  •  • 

02, iB  02,2B  •  •  ■  a2,jvB 


Z  =  A  (S>  B  = 


aA/,iB  ojvf,2B  •  •  •  om.jvB 
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Hadamard  Product:  A  and  B  must  both  he  M  x  N 

ai,2&l,2  •  •  •  Ol.JV&l.JV 

ry  A  <^2,1?>2,1  0‘2,2b2,2  •  •  •  0'2,Nb2,N 

Z  =  A  O  B  = 

_  0'M,lbM,l  0‘MflbM,2  *  •  •  0,M,NbM,N 

Vandermonde  Matrix:  given  an  M  x  N,  each  Xm  will  be  unique. 

1  xi  xl  •  •  •  ref 

1  X2  xl  •  •  •  ref 

■t  ^  /»2  ,»JV— 1 

1  Xm  Xm  •  •  •  Xm 
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Appendix  B.  Supplemental  Figures 


The  following  figures  show  the  full  adapted  angle-Doppler  pattern  for  a  target  at 
20°,  100  Hz,  and  each  scatterer  orientation. 
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Figure  56  Angle-Doppler  pattern  for  a  target  at  20°,  100  Hz,  with  no  scatterers. 
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Figure  57  Angle-Doppler  pattern  for  a  target  at  20°,  100  Hz,  with  scatterers  with 
(Ta  =  0.5  m^  placed  in  orientation  A.  The  number  of  scatterers  is  indi¬ 
cated  to  the  right  of  each  figure. 
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Figure  58  Angle-Doppler  pattern  for  a  target  at  20°,  100  Hz,  with  scatterers  with 
CTa  =  0.5  m^  placed  in  orientation  B.  The  number  of  scatterers  is  indi¬ 
cated  to  the  right  of  each  figure. 
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Figure  59  Angle-Doppler  pattern  for  a  target  at  20°,  100  Hz,  with  scatterers  with 
(Ts  =  0.5  m^  placed  in  orientation  C.  The  number  of  scatterers  is  indi¬ 
cated  to  the  right  of  each  figure. 
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Figure  60  Angle-Doppler  pattern  for  a  target  at  20°,  100  Hz,  with  scatterers  with 
CTs  =  0.5  placed  in  orientation  D.  The  number  of  scatterers  is  indi¬ 
cated  to  the  right  of  each  figure. 
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Figure  61  Angle-Doppler  pattern  for  a  target  at  20°,  100  Hz,  with  scatterers  with 
as  =  0.5  m^  placed  in  orientation  E.  The  number  of  scatterers  is  indi¬ 
cated  to  the  right  of  each  figure.  ^ 
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Figure  62  Angle-Doppler  pattern  for  a  target  at  20°,  100  Hz,  with  scatterers  with 
(Ts  =  0.5  m^  placed  in  orientation  F.  The  number  of  scatterers  is  indi¬ 
cated  to  the  right  of  each  figure. 
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